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ABSTRACT 

This  report  is  a  description  of  a  Fast  Fourier  Transform  (FFT) 
subroutine  written  in  assembly  language  for  the  ILLIAC  IV  computer.   The 
subroutine  uses  the  Cooley-Tukey  algorithm  for  performing  discrete  Fourier 
transforms.   The  parallel  nature  of  the  Cooley-Tukey  method  lends  itself 
very  well  to  a  highly  parallel  machine  like  ILLIAC  IV.   Timing  simulation 
results  have  shown  that  this  program  will  perform  Fourier  transforms  faster 
than  they  are  being  done  on  any  existing  computer  system. 
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Introduction 

The  Fourier  transform  is  one  of  the  most  important  mathematical 
aids  to  signal  processing.   The  basic  property  of  the  Fourier  transform  is 
its  ability  to  distinguish  waves  of  different  frequencies  which  have  been 
additively  combined.   The  result  of  a  Fourier  transform  on  a  complex  func- 
tion of  time  is  a  frequency  spectrum  of  that  function.   This  resulting  fre- 
quency spectrum  is  also  a  complex  function.   To  facilitate  processing  on 
digital  machines,  a  method  for  transforming  discrete  functions  of  time  into 
discrete  functions  of  frequency  has  been  derived.   The  time  function  must 
consist  of  evenly  spaced  samples  of  some  interval^  T.   A  function  of  N  samples 
could  be  represented  as  f(nT)  for  0  <  n  <  N  -  1.   The  spectrum  of  the  func- 
tion would  also  consist  of  discrete  values  evenly  spaced  in  frequency.   If 
the  spectrum  were  composed  of  N  points,  the  interval  of  frequency,  q,   would 
be  2tt/NT.   The  spectrum  F(koJ,  0  <  k  <  N  -  1  would  be  given  by 

N-l 
F(koJ  =  Z   f(nT)e"jnTkn  (l) 

n=o 

This  equation  defines  the  discrete  Fourier  transform  (DFT)  of  a 
sequence  of  N  samples.   It  follows  directly  that  the  DFT  is  linear  since 

DFT[f(nT)  +  g(nT)}  =  DFT{f(nT)}  +  DFT{g(nT)}  (2) 


and 


DFT(c(f(nT)))  =  c[ DFT{f (nT) } ]  (3) 

This  follows  quickly  from  the  definition. 

There  exists  an  inverse  DFT,  a  transform  which  maps  a  discrete 
Fourier  transform  back  into  the  sequence  from  which  it  was  composed.   It  is 
given  by 

f(nT)  .  |  lfo     FdB).*""  W 
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which  differs  from  the  DFT  equation  1,  only  by  a  scale  factor  and  by  the 

sign  of  the  exponential.   The  same  procedure  used  for  the  DFT  can  thus  be 

applied  to  a  frequency  spectrum  to  perform  the  inverse  Fourier  Transform  by 

simply  dividing  the  result  by  N. 

2 
The  above  method  for  computing  DFT  requires  N  complex  operations 

to  be  performed  in  order  to  transform  a  function  with  N  samples.   The  method 

of  Cooley  and  Tukey  [1]  eliminates  redundant  steps  in  that  method  resulting 

in  a  procedure  requiring  only  N  logp  N  complex  operations.   This  method  is 

highly  parallel  and  therefore  very  well  suited  for  use  on  ILLIAC  IV  [2]. 

Subroutine  Specifications 

The  ILLIAC  IV  fast  Fourier  subroutine  accepts  discrete  input  func- 
tions composed  of  N  complex  samples  where  N  is  an  even  power  of  two  between 
8  and  U096,  inclusive.   The  output  of  the  subroutine  is  a  discrete  frequency 
spectrum  composed  of  N  complex  samples  as  defined  by  equation  (l).   The  out- 
put is  stored  over  the  input  in  ascending  sequence  in  a  contiguous  block  of 
memory.   N  words  of  storage  are  used,  where  each  word  contains  the  real 
(outer)  and  imaginary  (inner)  portions  of  one  sample  in  the  32-bit  format. 
Input  to  the  subroutine  must  also  be  stored  in  this  manner;  therefore,  all 
values  have  2^-bit  significance  and  are  represented  in  floating  point. 
Several  functions  can  be  input  to  the  subroutine  at  one  time  as  long  as 
they  are  stored  contiguously  and  provided  that  they  all  are  the  same  length. 
If  the  transform  of  a  function  of  M  samples  is  desired  where  M  is  not  a 
power  of  two,  N-M  zero  samples  should  be  added  to  the  function  where  N  is 
the  smallest  power  of  two  greater  than  M.   The  result  can  be  input  to  the 
subroutine  as  a  function  of  N  samples.   The  output  will  be  a  frequency 
spectrum  equivalent  to  the  transform  of  the  original  function  of  M  samples 
with  a  smaller  frequency  spacing. 

Implementation  of  the  Cooley-Tukey  Procedure 

The  method  of  Cooley  and  Tukey  [1]  is  illustrated  in  Figure  1. 
Implementation  of  this  method  on  ILLIAC  IV  has  been  discussed  by  G.M.  Ackins 
[3]-   Ackins' methods  have  been  refined  and  generalized  for  use  in  this  sub- 
routine.  The  following  description  gives  the  details  of  the  final  imple- 
mentation. 
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Figure  1 


In  Figure  1,  f  represents  f(nT)  and  F  represents  F(npJ.   Also  W 

.2 
is  used  to  replace  the  constant  e  N.  At  each  node  in  the  diagram,  a  com- 
plex multiply  involving  the  lower  of  the  two  inputs  to  the  mode  and  the 
constant,  W,  raised  to  the  appropriate  power  is  performed.   The  result  is 
added  to  the  upper  of  the  two  inputs  giving  an  answer  which  is  used  as  an 
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input  to  the  next  iteration  until  the  transform  is  completed.   At  each  stage 
in  the  process  where  an  iteration  is  complete,  the  partial  results  are 
stored  into  the  initial  locations  so  that  the  data  is  transformed  in  place. 
Figure  \    illustrates  how  this  transform  would  be  accomplished  in  parallel 
on  an  eight  PE  ILLIAC  IV. 


Initial 

_data_ 

f 


First      Second  Third 

iteration   iteration  iteration 


Figure  2 
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In  Figure  2,  the  arrows  represent  routing  of  data  between  PE's 
while  the  nodes  represent  the  operations  which  take  place  within  the  PE's 
as  given  for  Figure  1.   f(t),  and  all  of  the  necessary  powers  of  W  are 
stored  in  PE  memory  for  use  in  the  calculation.   From  Figures  1  and  2,    it 
is  seen  that  the  input  data  has  been  "shuffled"  so  that  the  output  appears 
in  ascending  sequence.   This  shuffling  also  causes  the  routing  (arrow) 
pattern  to  be  very  regular.   The  algorithm  for  ordering  input  data  is  ex- 
plained in  a  later  section  of  this  report  which  describes  the  routine 
called  Scramble  which  performs  the  shuffling  in  ILLIAC  IV. 

The  method  described  above  and  illustrated  in  Figure  2  can  be 
generalized  to  transform  much  larger  functions.   Careful  examination  of 
Figure  1  reveals  that  after  two  iterations,  two  separate  four  point  trans- 
forms have  been  completed.   A  four  point  transform  is  depicted  in  Figure  3« 
It  should  be  noted  that  the  value  of  W  depends  on  the  size  of  the  function 
being  transformed. 
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Figure 

3 

The  eight  point  transform  can  now  be  represented  in  terms  of  the 
four  point  transform  as  depicted  in  Figure  h. 
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The  equality  of  these  two  representations  (Figure  1  and  Figure  k) 

is  important  to  the  structure  of  the  subroutine.   In  an  eight  point  transform, 

_.2  ^  _.  _.2 

W  =  e  B,  so  that  W  =  e   .   In  a  four  point  transform,  W  =  e  ^,  so  that 

W  -  e     .      Thus,  an  eight  point  transform  can  he  accomplished  by  performing 
two  four  point  transforms  and  iterating  once  more.   In  the  general  case 
any  N  point  transform,  where  N  is  a  power  of  two,  can  he  accomplished  by 
performing  two  N/2  point  transforms  and  combining  those  results  with  one 
final  iteration.   One  consequence  of  this  is  that  the  first  step  in  per- 
forming any  size  transform  is  to  perform  an  eight  point  transform  on  each 
group  of  eight  points.   The  entire  subroutine  is  built  up  from  that  basis. 

The  other  consequence  of  this  partitioning  of  the  procedure  is 
that  each  stage  is  completely  independent.   When  a  function  composed  of 
128  data  samples  is  to  be  transformed  on  ILLIAC  IV,  it  is  simply  stored  in 
two  rows  of  PE  memory.   One  row  of  PE  memory  consists  of  one  word  in  each 
PE.   The  first  6k   samples  are  stored  in  one  row  (after  proper  scrambling) 
and  the  second  half  are  stored  in  the  other  row.   This  is  referred  to  as 
"folding"  128  words  of  data  into  two  6k  word  rows  of  memory.   A  6k   point 
transform  can  now  easily  be  performed  on  the  first  half  of  the  data  since 
this  operation  is  independent  of  the  other  half.   This  6k   point  transform 
is  a  simple  extension  of  the  eight  point  transform  of  Figure  2,  using  6k 
PE's.   Another  independent  6k   point  transform  is  then  performed  on  the 
second  half  of  the  data. 

The  results  of  these  two  6k   point  transforms  can  be  combined  in 
one  simple  iteration  giving  the  final  128  point  transform.   Since  the  in- 
put data  was  folded  in  PE  memory,  no  routing  is  required  to  align  the  two 
halves  for  this  final  iteration  simplifying  it  even  further.   This  can  be 
illustrated  by  observing  that  the  final  iteration  of  each  transform  re- 
quires a  route  of  N/2  positions.   For  a  128  point  transform  this  would  be 
a  route  of  6k   which  is  equal  to  a  route  of  0  in  a  single  quadrant  ILLIAC 
IV.   Furthermore,  every  transform  larger  than  6k   requires  a  final  route 
which  is  a  multiple  of  6k   which  means  no  routing  is  actually  required. 
Since  there  is  a  smaller  proportion  of  time  spent  in  routing,  the  largest 
transform  is  actually  the  most  efficient  in  its  computation. 
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The  final  generalization  is  that  each  larger 'transform  above  6k 
is  simply  accomplished  by  folding  the  function  and  treating  the  halves  in- 
dependently until  the  final  iteration.   The  size  of  the  function  will  de- 
termine how  many  eight  point  transforms  are  done  as  the  first  step.   Also 
the  number  of  functions  to  be  transformed  at  a  time  can  be  taken  into 
account  through  proper  loop  control.   For  example,  two  20U8  point  transforms 
could  quite  easily  be  done  at  the  same  time.   They  would  each  be  scrambled 
separately  and  then  the  transform  would  proceed  as  if  the  two  functions  com- 
posed a  single  U096  point  function—only,  the  final  iteration  would  not  be 
performed. 

Scrambling 

The  shuffling  of  data  in  the  processing  elements  before  perform- 
ing the  transform  allows  the  computation  to  be  done  in  place,  that  is,  by 
writing  all  intermediate  results  over  the  original  data.   This  process  will 
be  referred  to  as  prescrambling  and  when  it  is  applied  the  transform  gen- 
erates results  which  are  sequentially  ordered  across  PE's.   The  original 
function  is  composed  of  discrete  samples  at  even  spaced  time  intervals. 
They  are  assigned  indices  0  through  N  in  order.   Each  index  can  be  rep- 
resented by  a  binary  number  with  log  N  digits,  KnK  . . .  K  ,  where 
n  =  log  N  -  1.   The  scrambling  process  requires  that  each  sample  having 

index  KAK.  . . .  K  be  interchanged  with  the  sample  having  index  K  K  ,  ...  K^. 
0  1     n  c  n  n-1      0 

Thus  by  a  simple  "bit  reversal"  of  each  binary  index  the  destination  of  each 
sample  in  the  scrambling  process  is  determined.   Note  that  in  Figure  1  a 
three  digit  bit  reversal  on  the  sequence  0,1,2,3A>5>6>7  yields  the  proper 
sequence  for  input  to  the  transform  procedure,  i.e.,  0,^,2,6,1,5,3,7-   Con- 
sider the  following  example  to  see  how  the  scrambling  is  accomplished  within 
ILLIAC  IV.   The  input  data  is  represented  by  indices  only. 
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PEo 

Initial 

function   0   1   2   3  k       5   6   7   8   9   10   11   12   13 

Bit  reversed 

indices    0   8  k     12   2  10   6  lU   1   9    5   13    3   11 

1st 

iteration  0   1   2   3   U   5   6   7   8   9   10   11   12   13 

Route  (7)     * *      t      T  ~*       t       * 


2nd 

iteration  0   8 

Route  (2) 

Properly 
ordered    0   8 


PE15 

lk       15 

7   15 

lk      15 


2  10   J+  12   6  li+   l   9    3   11    5   13   7   15 


k      2   2  10   6  lU   1   9    5   13    3   11   7   15 


This  process  can  be  generalized  for  scrambling  N  points  in  the 
following  manner.   The  odd  numbered  PE's  less  than  N/2  send  data  to  the 
corresponding  even  numbered  PE's  greater  than  w/2  and  receive  data  back 
from  them.   The  next  iteration  is  the  same  except  the  first  N/2  PE's  and 
the  second  N/2  PE's  perform  it  independently  and  handle  the  data  by  pairs. 
Note  that  only  two  routes  are  required  for  each  iteration.   By  this  method 
ILLIAC  IV  can  scramble  a  6k   point  data  set  in  three  iterations  requiring 
only  6  routes.   For  larger  data  sets  the  data  must  be  folded  in  storage  and 
a  different  scheme  is  used.   The  following  example  shows  how  a  sixteen  point 
data  set  would  be  scrambled  by  a  k   PE  ILLIAC  IV. 


PEo 

PE 

PE2 

PE3 

L0C 

0 
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3 

L0C+1 
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6 
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LOC+2 
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10 

11 

L0C+3 

12 

13 

lk 

15 

Proper  reordering 
determined  by  bit  reversal 
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13 
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11 
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15 
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The  underlined  indices  denote  the  positions  that  are  unaffected 
by  scrambling  since  data  in  these  positions  does  not  need  to  be  moved.   The 
location  of  these  stationary  positions  within  each  PE  is  given  by  perform- 
ing a  bit  reverse  operation  on  the  PE  numbers.   Each  PE  now  loads  a  register 
with  data  from  the  location  corresponding  to  the  stationary  location  in  the 
next  higher  numbered  PE. 


Register 
Memory 


L0C 

Lyic+i 

L0C+2 

LOC+3 


8. 

.       ^ 

11+ 
> 

3> 

PEo 

PE1 

PE2 

PE3 

0 

k 

12 

1 

5' 
9 
13 

2 

6 

10 

Ik" 

3J 

7 
11 

15 

Each  PE  routes  the  data  in  its  register  to  the  next  higher 
numbered  PE.   The  rightmost  PE  routes  back  to  the  first  PE. 


PEQ    PE     PE2    PE 


register 


Ik 


K7 — o — <cr 


PEQ    PE1    PE2    PE3 


register 


3 

8 

5 

Ik 

Now  the  data  in  each  register  is  exchanged  with  the  data  from  the 
location  in  the  PE  corresponding  to  the  stationary  location  in  the  next 
lower  numbered  PE. 
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PEQ    PE1    PE2    PE 
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Memory 
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L0C+2 
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9 
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15 

The  data  in  the  registers  is  new  routed  hack  one  PE  and  stored 
where  the  first  data  came  from. 


PE, 


PE1    PE2    PE 
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Memory 
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L$C+3 
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13 
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1J 
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Ik 
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15 
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The  transfers  accomplished  can  he  illustrated  as  shown  below: 


L0C 
L0C+1 
L0C+2 
L0C+3 


PEQ   PE]_   PE2 


PE, 


■11- 


The  transfers  to  complete  the  scrambling  would  be: 


PE 


PE. 


PE, 


PE, 
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12 
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Ik 

L0C+2 

1 

9 

5 

^11 

L0C+3 

.  3 

13-" 

"~T" 

15 

This  final  transfer  is  accomplished  by  only  one  route  since  each 
data  point  moves  exactly  two  PE's  away  (end  around).   Once  again  all  PE's 
are  active.   The  following  diagram  shows  how  this  final  step  is  accomplished 


PEQ  PE1 

PE2 

PE3 

L0C 
L0C+1 

0 

8. 

r2/ 

12 

*ir 

^LO 

6 

Ik 

L0C+2 

i 

9 

5^ 

~±±  —  - 

L0C+3 

3, 

*-13r^ 

7 

15 

1 

Thus,  after  two  iterations,  the  data  has  been  properly  ordered 
for  input  to  the  transform  routine. 
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12 

L0C+1 

2 

10 

6 

Ik 

L0C+2 

1 

9 

5 

13 

L0C+3 

3 

11 

7 

15 
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In  the  operations  described  above,  each  PE  is  always  performing 
some  tasks  and  the  data  is  never  moved  to  an  intermediate  memory  location 
so  that  the  scrambling  is  done  in  the  minimum  number  of  steps  and  with  100$ 
efficiency  in  the  use  of  PE's.   It  can  be  seen  that  neither  of  these  state- 
ments is  true  about  the  scheme  described  at  the  beginning  of  this  section 

or  the  scheme  described  by  Ackins  [3]-   Unfortunately  this  algorithm  works 

2 
only  on  data  sets  which  have  N  points,  where  N  is  the  number  of  PE's  in 

the  array. 

Since  ILLIAC  IV  has  6k   PE's  it  can  scramble  a  folded  data  set 
of  U096  points  100$  efficiently.   The  algorithm  has  been  extended  to  handle 
small  data  sets. 

Consider  the  following  example  of  scrambling  a  32  point  data  set 
in  an  8  PE  ILLIAC.   The  initial  state  would  be: 


L0C 
L0C+1 
L0C+2 
L0C+3 


PE, 


PE. 


PE, 


PE. 


PE, 


PE, 


PE, 


PE. 


0 

1 

2 

V 

k 

5 

6 

7  V 

> 

17 

10 

11 

*19 

V 

20 

21 

Ik 
^22 

15 
'23 

2k 

25 

26 

27 

28 

29 

30' 

31 

K_ 

\ 

Now  perform  the  transfers  indicated.   It  should  be  apparent  that 
these  operations  take  four  rather  than  two  routes  since  each  group  of  four 
PE's  cannot  be  treated  separately  in  an  end  around  fashion.   Also  these 
transfers  cannot  be  100$  efficient  since  more  than  one  routing  distance  is 
involved. 

The  result  of  the  transfers  would  be: 


L0C 
L0C+1 
L0C+2 
L0C+3 


PE, 


PE. 


PE, 


PE. 


PE, 


PEr 


PE, 


PE, 


0 

16 

^2 

2k 

k 

20 

^r6 

28 

8* 

-^T 

10 

26 

12  <4- 

■^22 

Ik 

30 

1 

17 

9 

J*1^ 

5 

21 

13 

JT23 

3 

25*- 

"*"ii 

27 

7 

29*- 

""15 

31 
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Now  only  the  above  indicated  transfers  are  required  to  complete 
the  scrambling  process.  Note  that  the  process  is  simply  two  applications 
of  the  previous  method  performed  in  parallel.   The  final  result  is: 


L0C 

L0C+1 
L0C+2 
L0C+3 


PE, 


PE, 


PE, 


PE, 


PE, 


PEr 


PE, 


PE. 


0 

16 

8 

2k 

h 

20 

12 

28 

2 

18 

10  ' 

26 

6 

22 

1U 

30 

1 

17 

9 

25 

5 

21 

13 

29 

3 

19 

11 

27 

7 

23 

15 

31 

A  generalized  method  for  scrambling  data  which  is  folded  in 
memory  can  now  be  given.   Data  is  first  scrambled  in  groups  which  are  com- 
posed of  as  many  FE's  as  there  are  rows  of  data  by  the  method  developed  for 

2 
"square"  data  sets  (referring  to  N  data  points  in  N  PE's).   Finally,  con- 
sider these  "squares"  as  units  and  scramble  their  position.   One  further 
example  will  make  this  last  step  more  clear.   Consider  16  data  points  in  an 
8  PE  array. 


First  iteration 


PEr 


PE. 


PE, 


PE, 


PE, 


PEr 


PE, 


PE. 


L0C 
L0C+1 


0 

1 

2 

-3 

h 

5 

6 

^1 

8*- 

9 

10*" 

^*D 

12  *- 

^13 

Ik** 

15 

The  heavy  lines  below  indicate  the  second  step  which  completes 
the  scrambling,  giving  the  proper  result. 


L0C 
L0C+1 


PE 


0 


PE. 


PE, 


PE, 


PE 


1+ 


PE, 


PE 


6 


PE, 


7 


0 

8 

2 

10 

k 

12 

6 

1-4 

1 

9 

3 

11 

5 

13 

7 

15 
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Result 


L0C 
L0C+1 


PE, 


PE. 


PE, 


PE, 


PE, 


PEr 


PE, 


PE, 


0 

8 

k 

12 

2 

10 

6 

... 
Ik 

1 

9 

5 

13 

3 

11- 

7 

15 

Although  the  present  implementation  allows  for  only  up  to  U096 
data  points,  larger  data  sets  could  be  handled  and  it  is  easily  seen  that 
scrambling  would  be  100%  efficient,  but  the  core  restrictions  on  the  size 
of  the  data  set  and  the  size  of  the  program  would  eventually  limit  the  pro- 
gram.  A  short  discussion  of  these  restrictions  is  included  under  future 
considerations  - 


Future  Considerations 

Various  extensions  to  this  subroutine  have  been  explored.   Some 
of  them  seem  impractical  while  others  are  quite  worthy  of  further  investi- 
gation and  perhaps  implementation.   There  also  has  been  a  suggestion 
that  a  separate  subroutine  be  written  to  handle  several  parallel  transforms, 
each  with  all  of  its  data  within  one  PE.   If  all  transforms  could  be  done  in 
groups  of  6k   this  would  appear  to  be  a  very  efficient  method;  however,  if 
groups  of  data  are  not  even  multiples  of  Gk,    some  PE's  will  be  completely 
idle  during  the  whole  process.   It  is  this  author's  opinion  that  such  a 
subroutine  would  be  very  straightforward  to  code  and  would  have  enough  value 
to  be  worth  programming  at  some  future  time. 

The  question  of  expanding  the  present  subroutine  to  accept  data 
sets  larger  than  U096  sample  is  largely  decided  by  available  core.   There 
does  not  seem  to  be  too  much  demand  for  greater  capability  at  this  time  and 
the  program  itself  already  consumes  approximately  2,500  words  of  memory, 
mostly  composed  of  stored  constants.   Each  increase  of  a  power  of  two  in 
the  size  of  data  sets  transformable  will  increase  the  size  of  the  subroutine 
by  almost  a  factor  of  two.   Such  increases  can  only  be  justified  by  a  large 
demand  for  that  capability.   If  and  when  expansion  is  warranted,  the  process 
would  be  straightforward  with  the  possible  restrictions  of  a  6kK   limit  on 
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the  size  of  the  data  set  so  that  it  could  be  completely  core  contained. 
Problems  with  formatting  and  accessing  disk  files  tend  to  be  formidable. 

It  has  been  shown  by  Robinson  [h]   that  two  real  data  sets  can  be 
handled  simultaneously  as  if  they  composed  one  complex  -data  set  of  the  same 
size.   It  was  suggested  that  capability  be  added  for  handling  real  data  in 
this  manner,  with  hopes  of  a  factor  of  two  increase  in  efficiency.   The 
algorithm  was  studied  and  appears  simple  enough  but  applying  it  to  a 
parallel  structure  like  ILLIAC  IV,  does  not  seem  feasible  at  this  time. 
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ILLIAC  IV  FAST  FOURIER  TRANSFORM 
USERS  MANUAL 

This  report  is  intended  to  be  a  description  of  how  to  use  the 
Fast  Fourier  Transform  subroutine  and  exactly  -what  it  can  do.   An  explanation 
of  the  Fourier  Transform  and  of  how  the  subroutine  works  on  ILLIAC  IV  is  given 
in  ILLIAC  IV  Report  #226  entitled  "A  Fast  Fourier  Transform  Subroutine  for 
ILLIAC  IV.  "  Also,  copies  of  the  assembly  language  listing  of  the  subroutine  are 
available  from  the  Center  for  Advanced  Computation  and  the  University  of  Illinois. 

The  Fast  Fourier  Transform  is  used  to  transform  a  discrete  time  series 
of  data  points  into  a  discrete  frequency  spectrum.   Both  input  and  output  data 
values  are  represented  as  complex  numbers  in  the  following  way:   each  number 
occupies  on  6U-bit  word  using  the  32-bit  floating  point  configuration  where  the 
real  part  of  each  number  occupies  the  inner  32  bits.  When  the  input  series  con- 
sists of  only  real  numbers  the  numbers  should  be  stored  only  in  the  outer  por- 
tion of  each  word  and  the  inner  portions  should  be  set  to  zero.   In  any  case  the 
output  will  always  consist  of  complex  numbers  and  will  be  stored  in  the  same 
locations  as  the  input,  thus  "overwriting"  the  input. 

Due  to  the  structure  of  the  Fast  Fourier  Transform  algorithm  all  input 
sequences  must  contain  a  number  of  data  points  which  is  an  even  power  of  two. 
The  sizes  of  data  sets  which  can  be  transformed  by  the  ILLIAC  IV  subroutine  are: 
8,  16,  32,  6k,    128,  256,  512,  102^,  20^8,  ^096.   Other  sizes  of  data  sets  can 
be  transformed  by  "padding"  them  with  zeros  on  the  end  until  their  size  is  in  an 
even  power  of  two  between  8  and  U096  inclusive.   The  result  of  such  a  transform 
will  be  correct  with  only  the  frequency  spacing  of  the  result  being  effected. 

This  subroutine  uses  the  standard  subroutine  linking  conventions  as 
described  in  "The  ILLIAC  IV  Software  Reference  Manual. "  When  referenced  from  a 
GLYPNIR  program,  a  standard  subroutine  call  statement  can  be  used.  When  refer- 
enced from  an  assembly  language  program,  the  standard  call  macro  can  be  used. 
This  macro  and  an  example  of  its  use  is  listed  in  Appendix  A.  The  Fast  Fourier 
Transform  subroutine  requires  three  parameters  named  size,  number,  and  pointer. 
These  parameters  have  the  following  significance. 

SIZE 

Size  gives  the  number  of  data  points  present  in  each  input  time 
series.   Size  must  be  an  integer  and  should  be  an  even  power  of  two  between  8 
and  U096  inclusive.   If  it  is  not  an  even  power  of  two,  the  next  greatest  power 
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of  two  will  be  assumed  and  will  determine  the  size  of  the  transform  and  the 
number  of  results  returned.   Note  that  the  data  must  still  be  "padded"  with 
zeros  for  the  result  to  be  meaningful.   If  the  size  specified  is  zero  or  greater 
than  4096,  an  error  condition  will  be  indicated  and  no  transform  will  be  per- 
formed.  An  error  is  indicated  by  filling  ACAR2  with  all  ones. 

NUMBER 

Number  must  also  be  an  integer  specifying  the  number  of  transforms 
to  be  performed  during  this  call  to  the  subroutine.   Any  number  of  transforms 
can  be  performed  as  long  as  they  are  all  the  same  size  and  as  long  as  all  of  the 
input  data  ia  available  in  one  continuous  block  of  memory.   If  number  is  set  to 
zero,  no  transform  will  be  performed  and  ACAR2  will  be  set  to  all  ones. 

POINTER 

Pointer  is  a  2^-bit  CU  address  giving  the  location  of  the  first  word 
of  a  continuous  block  of  memory  which  contains  all  of  the  input  data  for  this 
call  of  the  transform  subroutine.   Pointer  will  also  address  the  first  word  of 
the  resulting  frequency  spectrum(s)  when  the  subroutine  is  exited.   This  pointer 
is  automatically  set  up  by  the  calling  macro  which  is  exhibited  in  Appendix  A. 

The  actual  subroutine  is  contained  in  a  file  called  FASTFOURIER/ 
TRANSFORM  which  is  available  on  a  system  tape  of  the  subroutine  library.   This 
file  must  be  properly  linked  to  the  calling  program.   The  subroutine  is  called 
by  using  the  entry  identifier  FFT. 

INVERSE  TRANSFORM 

For  performing  inverse  transforms  the  only  change  required  is  to  call 
by  using  the  entry  identifier  FFTINV. 
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APPENDIX  A 

Standard  call  macro  from  "The  ILLIAC  IV  Software  Reference  Manual" 

DEFINE  CALL  &NAME(&PARAMETERS)= 

&IF  &SIGN(&MFIELD  (&NAME) ) 

&THEN  EXTERNAL  &NAME;  &FI 

BEGIN  BLOCK 

BEGIN  USE  63 

LIST:   DATA  &PARAMETERS; 

END; 

CLC(2); 

SLIT(2)   LIST; 

END; 
clc(3); 

SLIT(3)   &NAME; 
EXCHL(3)   $ICR;  ##; 

Sample  calling  sequence  to  execute  two  I38  point  transforms. 

CALL  FFT(HOWBIG,HOWMANY, WHERE); 


WHERE:   DATA       [128  words  of  data]; 
HOWBIG:   DATA      6k; 
HOWMANY:   DATA     2; 


■19- 


APPENDIX  B 
Demonstration  Run 

This  appendix  compares  a  Fourier  Transform  performed  on  a  Univac 
1108  and  the  same  transform  performed  on  the  ILLIAC  IV  Simulator.   The  infor- 
mation is  presented  in  the  form  of  actual  computer  listing  including  the 
following:   l)   A  listing  of  the  256  real  input  data  points.   2)  The  results  of 
the  Univac  run  (real  part  first  followed  by  imaginary  part).   3)  The  main 
program  used  to  call  the  FFT  routine  for  the  ILLIAC  IV  simulation  run. 
k)      The  results  of  the  simulation  run.   Differences  in  the  two  results  are  due 
primarily  to  the  difference  of  the  word  size  used  in  each  case. 
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Input  Data 

o. 

0. 

-5b. 

-24. 

4822. 

-13007. 

-9326. 

5835. 

171.67. 

1  2  D  4  7  , 

3. 

-696. 

-1676. 

77  34. 

-2215. 

-7602. 

1 0 21 1 . 

lo5. 

-17335. 

9888. 

-12968. 

2744. 

14245. 

10537. 

11^*1  3. 

l'WW, 

4049. 

<069. 

)  7389. 

6574. 

18416. 

507. 

|C625. 

2  364. 

-3213. 

-36  09. 

-6436. 

-11208. 

3291. 

7014. 

-7/'<0. 

-3226. 

-16879. 

1217  3. 

-684  8. 

-8584V 

7979. 

13897. 

12/3. 

- 1  5 1 B . 

1 1 723. 

3/09. 

768. 

-7109. 

3905. 

-17075. 

17524. 

-6860. 

8704. 

1967. 

-4606. 

-11202. 

16318. 

3432. 

4807. 

-19886. 

-6603. 

-11205. 

-11151. 

762. 

7412. 

951. 

31 1>7. 

32  70. 

11954. 

3  0  0  3  . 

-1335. 

-14578. 

17089. 

"-8457 

-6505. 

-12545. 

4537. 

5618. 

-5144. 

-13695. 

1017  0. 

-3683. 

-8/o2. 

-924. 

-2680. 

4/16. 

-4320. 

-5671. 

-11186. 

12147. 

-6bl  . 

1 1049. 

-12576. 

3387. 

1053. 

-12898. 

5476. 

1719. 

-5  6  92. 

-9780. 

611b. 

4889. 

-15757." 

""6186.  " 

-1947. 

"  -1902. 

9280. 

40  03. 

-13438. 

1  34  66. 

-14  0  7. 

-5415. 

2584. 

-13501. 

4027. 

-41 76." 

223b. 

-10937. 

85  7  3. 

"-13987." 

"3345. 

-1129. 

-530?. 

-9614. 

-64  15. 

-1  3410. 

1344  7. 

1673. 

-10575. 

-3613. 

]o5. 

-191. 

499. 

1514. 

-5063." 

"19279. 

"8  32. 

-5213. 

12544. 

49  /6. 

-17782. 

-H99fa. 

10533. 

-5940. 

-2264. 

7824. 

6105. 

-16655. " 

9  0  24/ 

3  7666. 

3006. 

-7317. " 

~  -8  118." 

3976. 

"5297. 

-109  75. 

10  253. 

16419. 

2224. 

12068. 

392. 

1382. 

~  - 3  1  2  . 

-238. 

1088b. 

1551. 

2753." 

6381  . 

"8903. 

-393. 

14*  ;23. 

11118. 

-9271. 

-7248. 

1997. 

-1026. 

809. 

4387. 

-3143. 

6325. 

-4  87  7. 

195. 

-496. 

17  32. 

16910. 

-3611. 

-5082. 

-3918. 

-5548. 

-2748. 

5597. 

3206. 

-8721. 

-7531. 

ill?. 

15922. 

293b.  " 

-  4  695." 

-2051. 

18  043. 

"""2013. 

"  -306. 

-5539. 

-901 7, 

-18840. 

10401. 

3189. 

-18287. 

-3289. 

-4301. 

-2925. 

-3157. 

10105. 

"3154. 

-7270. 

-.1  1441. 

"""2041  7" 

"-37467 

3139, 

-8121. 

-12835. 

9925. 

2178. 

-6024. 

-8531. 

4158. 

j-486]  . 

-8923. 

18110. 

13781. 

"3  060/." 

51297" 

157877" 

34287 

n-1069. 

-2 Ob 8. 
2816. 

-12438. 
1816/. 

1 2253. 

1954. 

-7699. 
-279  7. 

-10595. 
-6039. 

6591.. 

-1159. 

-15/5. 

17566. 

2517. 
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Real   Univac   1108   Result 


.85580000*04 

-.23935650*05 

.  10475777*06 

.40882889*05 

-.79964560*04 

.  16432617*06 

-.H81901?6*05 

. 12018498*06 

-.928207(12*05 

-.213009^2*05 

."3661  06P6.05 

.601555^6*05 

.  10175  7^4*06 

.23466077*05 

. 126766O5.06 

-.50596355*05 

.11267843*06 

.576140)0*05 

.28340147*05 

.930267^9*05 

-.12314756*06 

.21500096*06 

-.195869?4*06 

.39295759*05 

.17858835*05 

.123053h9*06 

.30633514+05 
-.79698125+05 
-.89308050+05 
-.11257891+05 

.11725652+05 

.19196831*06 
-.88314708*05 

.12051347+06 
-.29508836+05 

.54568208*05 
.50177875*05 

-.81431963+04 
-.65641566+05 

.62304680*05 

-.75575159*05 

.36708912*04 

.86965552*05 

.95429297*04 
-.94544403+05 

-. 17295822+06 
.17787883+06 

.54710958*05 
-.72483458*05 

-.28900264+06 
-.42655679+05 

-.66131058*05 
.11841725*06 

-.40119934+05 
.13246392+06 

.41877700*05 
-.10685043*06 

.35511073+05 
-.34118988+05 

.10627192*06 
.25315634*05 

-.61932731+05 
-.93928612+05 

.28277265*05 
-.86787733*05 

-.1  13600  16*06 
.63162532*05 
.62923921*05 

-.19840547*04 

-.44793332*05 

-.35014947*05 

. 19538353*06 

.6791216R+04 

-.86787713*05 

.28277265*05 

-"738144223*05 
-.84544034*05 

.67912168*04 
.19538353*06 

-.22494207*06 
.59866334+05 

-.35014947*05 
-.44793332*05 

.32308614*04 
.19907523*06 

-. 19840547*04 
.62923921*05 

.59887308*05 
-.25588855*05 

.63162532*05 
-.11360036*06 

-.24503150+05 
.75343836+05 

.12305369*06 
.17858835*05 

.25315634+05 
.106271^2*06 

-.11576040+06 
-.56760156+05 

.39295759*05 
-.19586924*06 

-.10685043*06 
.41877700*05 

-.49808327+05 
-.12489943+06 

.21500096*06 
-.12314756*06 

.11841725*06 
-.66131058*05 
-.72483458*05 

.54710958*05 

-.35125625+05 
.55541602+05 

.93026729*05 
.28340147*05 

.95541438*05 
-.10801400*06 

.57614010*05 
.11267843*06 

.86965552*05 
.36708912*04 

.18504945*05 
-.11136791+06 

-.50596355*05 
.12676695*06 

-.75575159*05 
.62304680*05 
.50177875*05 
.545682(18*05 

-.88314708+05 
.19196831+06 

-.79964560+04 
.40882889+05 
.10475777+06 

-.23935650+05 

-.12558822+00 
-.12877611+06 

.23466077*05 
.10175734*06 
.60155586*05 
.36610626*05 
-.21300992*05 
-.92820702*05 

.59602295+02 
.62851 124+05 

-.14447113+05 

.68437002+04 

-.12452422+06 

-.18619060+04 

.12018498*06 
-.88190126*05 

.12423822+06 
-.43274833+05 

.12878585+06 

•.32208377*05 

•.31424301*05 

-.14610218*05 

•.11268548+06 

.46726594+05 

.10651084*06 

-.35035923*04 

-.1 1844547*06 

•.22496512*05 

.23490922*05 

.24522155+05 

-.13490052+06 

-.70621265*05 

-.87760739*05 

.63107815*05 

_L18619414*05_ 

■.fB944333*06 

^87388329*05 

.48S61633*05 

-.35586001*05_ 

-.20365177*05 

.14108025*06 

.12"12353b"*06~" 


^.42466545*05 
.11306600*06" 
^.42466545*05 
.12123530*06 
.14108025*06 
•.20365177*05 
L.255R6001*05 
.48561633*05 
^87388329*05 
•."18944333*06 
^18619414*05 
.63107815*05 
^87760739^05 
•.70621265^05" 
••13490052*06 
".24522  155*05 
^23490922*05 
•.22496512*05 
..11844547*06 
••35035923*04 
.10651084*06 
.46726594*05 
•.1J268548*L06_ 
••14610218*05 
L.J  1424  30  1_*05 
■.32208377*05 
.12878585*06 


-.43274833*05 

.12423822*06 

-.18619060*04 

-.12452422*06 

.68437002*04 

-.144471  1  3*05 

.62851 124*05 

.59602295*02 

-.12877611*06 

-.  12558822*06 

■.11136791*06 

.18504945*05 

-.10801400*06 

.95541433*05 

.55541602*05 

-.35125625*05 

-. 12489943*06 

-.4980832"7*05 

-.56760156*05 

-.11576040*06 

.75343836*05 

-.24503150*05 

-.25588855*05 

".S"988730S*05 

.19907523*06 

.32308614*04 

.5986633^*05 

-.22494207*06 

-.84544034*05 

-.38144223*05 

-.93928612*05 

-.61932731*05 

^.341  18988*05 

.35511073*05 

.13246392*06 

r.46119934*05 

-.^♦2655679*05 

r. 28900264*06 

.17787883+06 

:.  17295822*06 

-.94544403*05 

.95429297*04 

-.65641566*05 

-.61431963*04 

-.29508836*05 

.12051347*06 

.11725652*05 

-.11257891*05 

-.89308050*05 

-.79698125*05 

.30633514*05 


.16432637*06 
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Imaginary  Univac    1108   Result 


.00000000 
.  12539  7-+9  +  05 

-,71297144+OS 

.  18855b8b+06 

.17435297*06 

.35799352*05 

-.626712  31*05 

-. 1 1224Qb9*06 

-.92065512*04 

-.83311300*05 

-.21239087*05 

-.78648c>52*04 

-.11474280*06 

.10789059*06 

-.57622153*05 

.18087399*05 

-.17384096*06 

-.97109678*05 

-.38454852*05 

-.48075970+05 

_.30471259*05 

.11034351*06 

.56991673*05 

-.731532/7*05 

-.14147976*06 

-.20342460*06 

.34484995*05 

-.75169292+05 

.73496373*05 
~.54843035+05~ 
_ .16088750*05 

.10316574*06 

.22981703*06 
-.20013107*06 
-.15650693*06 
-.50385053+05 
-.88054139*05 

.44369014+05 
_ .33982734+05 
-.10571027+06 
-.12952341*06 

.13493994+04 
_ .72885594+05 

.88596548*05 
-.46812337*05 
-. 14840  741 +06 

.12412578+05 
-.21966100*05 
-.299284)9*06 

.15715059*06 
-.94287190+05 

.32348535*03 


-.32348535+03 

.942H71'<0*05 

-.  15715054+06 

.2992H4 19+06 

.21966100*05 

-.  12412578+05 

. 1484074] +06 

.46812337+05 

-.8859b54M+05 

-.72885544+05 

-. 1 3493994+04 

.12952341 +06 

.10571027*06 

-.33982734+05 

-.4^369019+05 

.880541 19+05 

.50385053+05 

.15650693*06 

.20013107+06 

-.22981703+06 

-.10316574+06 

-.16088750+05 

-.54843035*05 

-.73496373+05 

.751892'-'2  +  05 

-.34484995+05 

.20342460+06 

.14147976+06 

.73153277+05 
-.56991673+05 

-.11034351+06 

-.30471259+05 

.48075970+05 

~  .38454 852+ 05 

.97109678+05 

. 17384046+06 

-.  18087399*05 

.57622153+05 

_-, 10789059+06 

. 1  1474200  +  06 

^.78648652  +  04 

.21239087+05 

.83311300+05 

.92065512+04 

.11224059+06 

.626  71231+05" 

-.35799352+05 

"-.17435297  +  06 

-.188556«6  +  06^ 

.71297144+05 

-.12539749+05 


-.15282431+06 
.34464606+05 

-.96626729+05 
.47544088+05 

-.10578316+06 

-.48ifl9366  +  05 
.649649/7+05 

-. 10760260+06 
.44764577+05 
.95072141*05 
.19846301*06 
.93203335*05 
.23333460+05 

-.15458620+05 
.35402116+05 
. 15852791+06 
.23991060+05 

-.79624333+05 

-.25320520+05 

-.14226242  +  06" 
^36430140+05 
.11887141+06 

-.74552653+05 
.188907)8+06 

-.92939891+05 
.12785482+06 
.33364525+05 
.19050124+05 

. 12846535+06 
.44  664971  +04" 
-.54591791+05 
-.74105499+05 
.48101416+04 
"  ."19423950  +  05 
.63503576+05 
.47605405+04" 
-_.  17098704+05 
-.28431246 ♦05 ~ 
-.23472000+06 
-.  10436177+06 
_-.44790«95  +  05_ 
."72368868  +  05 
-.94768170+05 
-.45146798+05 
.21326760+05 
-.36626935+05 
.32980707+05 
.76239434+05 
.20518874+05 
-.60283592+05" 
_-, .79  2  34  054  ♦05. 


-.27383049+05 

-.41893763+05 

.51  172803  +  05 

.49064532+05 

.1  1941343  +  05 

"."4696  7453  +  05 

-.17344318+06 

-.13977412+06 

.47270071 +05 

.24357188+06 

.31240915+05 

.14410096+06 

-.322141 13+05 

-.24044316+06 

-.31654141+06 

-.22730289+06 

-.13077045+06 

".  10901463  +  06 

-.47217832+05 

.32640422+05 

-.12036483+06 

.405873 19< 05 

-.11836559+06 

.12702667+06 

. 10989458+06 

-.00000000 

•. 10989458+06 

■. 12702667  +  06 

.11836559+06 

-"."40587319  +  05 

.12036483+06 

-.32640422+05 

.47217832+05 

-.10901463+06" 

^13077045+06 

.22730289+06 

.31654141+06 

.24644316+06 

.322141 13+05 

-.  14410096  +  06" 

-.31240915+05 

-.24357188+06 

-.47270071+05 

.  13977412  +  06 

_.  17344318  +  06 

-.46967453+05 

■^.11941393  +  05 

-.49069532+05 

-^51172803+05 

.41893783+05 

.._J_7383049_+05 


.79234054*05 

.60283592*05 

-.20518874+05 

-.76239434+05 

■^32980707  +  05 

.36626435+05 

-.21326760+05 

.45146798+05 

.947681 70+05 

-.72368868+05 

.44790895+05 

.104361 77*06 

_.23472000  +  06 

.28431246+05 

.  17098704  +  05 

-.47605405+04 

-.63503576+05 

■"."  19423950*05 

-.48101416+04 

.74105499+05 

.54591791*05 

-.44664971+04 

-. 12846535*06 

-.19050124+05 

-.33364525+05 

-.12785482+06 

.92939891+05 

-.18890718+06 

.74552653+05 
-.11887141+06 
-.36430140+05 

.  14226242  +  06 

.25320520+05 

,79624333+05 
-.239Q1060+05 
-.  15852791  +  06 
-.354021 16+05 

.15458620+05 
-.23333460+05 
-.93203335+05 
-. 19846301+06 
-.95072141+05 
-.44764577+05 

.10760260+06 
-.64964977+05 

.48219366+05 

.10578316+06 
-.47544088+05 

.96626729+05 
-.34464606+05 

..15282431  +  06 
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Main  Program  of  ILLIAC  IV  Simulation 

(Data  stored  in  external  file  referenced  by  ABC) 


BEGIN 


EXTERNAL 

FFT,ABC; 

SIZE: 

DATA 

256; 

NU: 

DATA 

i; 

MAIN: : 

SETE 

E.OR.-E; 

SETE1 

E.OR.-E; 

SLIT(2) 

ARGLIST; 

CLC(3) 

SLIT (3) 

FFT; 

EXCHL(3) 

$ICR; 

lit(o) 

I.ABC+6OO  .ABC 

DISPLAYR 

$C0.20:8; 

HALT; 

ARGLIST: 

DATA 

SIZE, NU, ABC; 

■2k- 


Illiac    IV   Simulation  Results 

(Outer   =  Real,    Inner   =   Imaginary) 
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APPENDIX  C 
Timing  Sunmary  t 


SIZE  OF 
DATA  SET 

1+096 

201+8 

102*+ 

512 

256 

128 

6k 

2  x  32 

1+  x  16 

8x8 

Overhead 


TIME  FOR 
SCRAMBLING 

75«0  ju  sec 

kS.h  u  sec 

53*1  u  sec 

23.2  \i  sec 

17.^  /u  sec 

10.5  u  sec 

I+.65  ju  sec 

2.80  in  sec 

2.1+0  u  sec 

1.75  jLi  sec 


TIME  FOR 
'FFT 

2233.6  (U  sec 

IO38. k  u  sec 

1+80.0  ju  sec 

220.1+  ju  sec 

100.1+  u  Sec 

1+5-35  u  sec 

20.6  ju  sec 

15.65  ju  sec 

11.10  in  sec 

6.75  U-  sec 


TOTAL 
2308.6  n  sec 
IO96.8  ju  sec 
533»1  u  sec 
2I+3.I+  ju  sec 
II7.8  ju  sec 
55.85  u  sec 
25.25  /i  sec 
18.1+5  u  sec 
13.50  ju  sec 
8.50  ju  sec|i 
13       i±  sec 


3+     x  ^ 

2301+  ju  sec 

1056  ju  sec 

1+80  a  sec 

216  ju  sec 

96  u  sec 

1+2  ^  sec 

18  ju  sec 

15  ju  sec 

12  ju  sec 

9  u  sec 


This  formula  gives  a  rough  estimate  of  the  time  required  to  do  the  FFT 
portion  where  each  complex  multiply  operation  takes  3  U  sec. 

t     All  times  base  on  50  n  sec  clock. 
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APPENDIX  D 

SUBROUTINE  LISTING 
LISTED     nM     0H/1J/71        AM     OV  »  14 

<E'-IN     CljMMtNT 

THE  FAST  rnijHIFW  Ihansfdhm  SUBROUTINE  IS  COMPOSED  OF  FOUR 
SL  ARATE  LUAD  MiiOULtS.    THE  FIRST  IS  THE  MAIM  PROGRAM   HlCH  HANDLES 
TMr  SUBROUTINE  |  IN  '  h  f  i  F  AN|l  SFIS  UP  ThE  C'l  REGISTERS.   ^F -T  COMES 
TMt  SCRAMBLING  ROUTINE   HlCH  REORDERS  THF  lM»Uf  nATrt  SO  TmAT  IT 
rAm  .;E  TMANSFH  1T|).   THEN  CUflES  THE  TmANSfJRM  routine  which 
CU-  PLETLS  THF  TRANSFORM,  PERFORMS  AM  IN'FRSF  IF  DESEED,  AND 
rEtURNS  rUNTROL  TU  T  HE  CALLING  UROGRAM.   FINALLY  ALL  UF  ThE 
C,U  STANTS  USED  ri|K  ThE  TRANrflHM  ARE  COMPILED  IN  A  SE»ARATF  HATA 
Ell  E. 

THREE  HARA»«FIEkS  AHE  RF-'UIPEO  HY  THIS  PRUliRAM,   THEY  ARE 
Li?TFU  HELO*  IN  THE  riRUF'R  Im  ''HTCH  THEY  MUST  HE  SUPPLIED. 

1)  ST^F  -  THF  HiimHEH  UF  HA  I  A  POINTS  PER  DATA  SET. 

?)     NUMUFR  "  THF  NUM-.F.H  OF  DATA  SETS  T<J  HF  TRANSFORMF D  . 

3)  POINTER  -  ThF  CU(?a  hTT)  AoorFSS  OF  TmE  FIRST  LOrATTnN  Or  DATA, 

A  USERS  MAMijAL  rUR  THlS  ill  ROUTINE  HAS  »EtN  PUtL.I$HEO  AMI)  IS 
TItLFU  ?THE  ILLTAC  IV  FAST  FOURJFH  TRhNSFURM  SU«HUUT  I NE  USEPS 
MA-  UAL"  HY  JAMF*;  E  STEVENS   It. 

ThF  main  pur  1 1  on  jy  ThE  pwoc.ram  FOLLOWS  ANn  SETS  Up  ^u 
Rt^i.sTERS  VITH  CONIRIIL  INFORMATION  As  SHO*N  HE  LH  j  • 
LUrAL  DATA  riUEFrR  CL'>M  OR  AO  ) 

DUf PATTERN) 

In   the   scramble      sErrhH   or    Tmf   p^nr.RA  1  This   REgIstEr   ts 

UbrU    F'lR    SAy/ING     MODE     PATTERNS 
0 1  r  F  L  A  G 1 

this  register  is  dwioe.j  into  three  i-ifl'>s.  tme  ftpst.  hit 
nu  >er  is.  is  used  ro  indicate  -hETher  an  inverse  transform  has 

FJtFN  RE'JUESTFi),   IF  -'IT  15  IS  SET  TU  ONE  AN  INVERSE  TRANSFORM  wTl.L 
mE  PFKFORMEO,   THE  SFCOMM  FTELO  OSES  THF  NEXT  in    HITS  A m i>  Is  SIMPL* 
A  rUPY  UF  THF  FjwST  PARAMETER  Sl'E,  WHTCH  GIVES  THF  NUMHFR  OF  OAT« 
PllTNTS  IN  EACH  TNHiJI  DATA  SET.   ThE  VARIABLE.  SI/E,  MUST  «E  AN  EVE' 
PU  Ep  f)F  T«U,    TF  THE  PARAMETER,  SI/E.  PASSED  TCI  THF  SUBROUTINE  IS 
NUt  A  PH.iER  HE  T"iU  IT  ^  I  L  L  h'  ASSUMED  THAI  TmE  NFXT  LARgEsT  PUwEP 
OK  JO  IS  UESIRFD  A-i)  THE  PROHf;n  VALllE  IS  THEN  GFNEpATEd.   ThE 
FI'AL  2H     *ll    FKLL)  UF  THE  r'L'A"  REGISTER  TS  LOADED  <"JTH  THF  PF 
AU'XFSS  OF  THE  rlRSF  LJCATiUn  IK  T  AV     DATA  SET  Tl)  HE  TRANSFORMED. 
THfS  AnnRESS  IS  E^OAL  M  THE  THTRI)  paRAMFIEw.  PniNTFR.  OlVlDEO  BY 
b*    SO  IT  IS  A  PF  ADDRESS. 

LATER*  IN  THE  TRANSFORM  SECTION  U^  THIS  PROGRAM,  ThE. 
INrllRMATIUN  In  REGISTER  DlULATo  IS  LOAOFU  INTU  ACAwO  AkjD  RFfERED 
TU  ny  THF  NAME  <U/F. 

D*(KFT(IRN) 

THIS    KE^ISTFH     IS    USED    T  '    SA>'E    ThE    RETURN    ADDRESS    EniiND    TN 
AC-R1     ^HFN     THE     SIMRuUTIvE     iS    ENTERED,        THE     ACTUAL    RFTU»M     ALsn    USFS 
AC/>Hl,        AN    ERROQ    KETiIKN     15     FLAG^FU    '^Y     FILLING     ACAR?     WITH    ALL     UNES. 

DJ     AND     D4 

SAVE  AREA  niH  ACARD  AND  AC.ARl. 

AC  f  RO( ST  71 .  INDFV  ) 

IN  THE  MAIN  SECTION  UE  THE  SOmRdUTInE  THIS  REGISTER  IS  L"ADEO 
WItH  ThE   VALUE  ST^E  «.  THF  FIpSI  PARAMFTER).   AFTER  T  hF  ST/F  TS  USED 
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TU    1AKE     a     CH  If'!  t  in  Jul      .-UA-JCh     I  -ITS    HtfilSTTN     IS    USED    Tfl    ^HLD    TNOExIN'. 
I»Nt  UHMATIDN    i)URT''<i     Ink.     5 C rt A  -i i- L I     •',    RPilCESS*       AS    NOTEn    A*n»/F»     FOP 
Tur     TRANSFORM     St  C  I  I    N      '►      Till',     pP.MiRAM     TuTS     RFGlsTtR     TS     |I)AUF'>     FROM 
Rt'.l5TF.R    i)1(FLAi;)     Aid    Pt  MIKF.i     !   i    AmATN    PY     TMF     NAME    SI*F, 

AC  «  P1  (  Mi). it  ) 

ThTS     Hrr,ISTn<     I  S     NIT     MlLFn     vITh     AMY     INITIAL     »/ALUF     AnU     TS     USE' 
THnOljdHflilT      TmP      ^iiHridijT  I-iF     f    W     S'WlNfi     MOI)F     PATTERNS     FHP     FNAU|.IN(,     pE  5 

ac  » w?( l hop  J 

TMTS     KF'ilSTf.'K     IS     '"  IT     MLLFn     /'I  T  H     AMY     I  N  T  T  I  A  |_     VALUF     Am()     TS     USE 
TH'j  JOl>Hni|I     THE     S'l^.lijTI    ^      Mi'<     U  riP     CONTROL     iNriiUMATTliN. 


A  C  '  k  i  (  c  n  1 1  n  T  ) 

TMTS    PF    .TSTf^    CiNTaTNb    4     PufNTFp     M     THE     RU*    OF    fl  A  T  4     ML  I  mG 
CU.  M  Ut'RFH*     AS       FLL    AS    i\    PUHIEP    Ti)    THE    LAST    PU*    C  j  f-     OATa     (A    L I  ^  I T  )  ." 
IT     Is    IJSFD    TO     I'ii)t<     pL     "■    1UPV     A'.»    M    OF.  Trp  M  T  NE     riHEN    THE     ToAnjsFuPM 
TS    COMPLFTt.        TmTS    RFi,  IsU'P.     ]5    sf  T     up    py    THE    mATn    SFCTlnN    Or     ThF 
SU  'NlTUTT'lE     AS    FniLH.SS     (1)     Tut      T'CmFmENT     FTFlO    IS    SFT     Tn    nNF  . 
(2       Tnf     L1"11    r'rL'i    TS    s1"  T    Tn    TJ    AonWEsS    Of    T  HF    FlwSl     ph..    nr    -iATA 
PL    S    ThE     TUfAL      l'10-H-.H    t'F     P'1  <S    S'IPPLTFO    (sTpF     Sl^F    fir    T^f    nATA    SETS 


1  I   ,ts     THF     NLT1  •*•  w     UF     1)6  t 
(LAST    2 "     HITS)     t  <;    s  L  I     T 


.SET-.     :iT./IOE|)     MY     64).     ID     Tpr     AnOoESS     FJFL'> 
Trif      i'F      APDPFSS     Of     TMF     FTPST     HOW    0?     nATA: 


t  X  T  F  R  N  A 1 

-EPt  : 

fill. 

1  '1 . 

FFTr  FN  K  Yi  »             SLT  T  (  n  1 

1) : 

SKIP 

.  1  s 

FFTINV/UNTKYI  :     Sl.TT(nj 

l  ; 

STL(O) 

h"i; 

STL(  1  ) 

♦•i;'J; 

STL(  n 

'•-  )  ?.  1 

LUAOC2) 

h  C  J : 

LOAif  *  ) 

*C  1 : 

Ai.iT(.n 

-1  ; 

CSHI.  (  i  ) 

1  ' 

LFA.XK  J) 

• 

CTSMF I  3  j 

'■>  !">  t  F  m  P 

c  p  p  r  n 

S4!  J 

SLITC1  ) 

*  H  5 

C$IIH(  1   ) 

%C  ^  : 

SL.IT(J) 

=  1  ; 

CSHI_(-<) 

(  l  >; 

CSHL(O) 

py : 

CA')O(U) 

«,  c  ■«  •• 

ALTT(^) 

1  ; 

LUAi)(  <>  ) 

s  c  .* ; 

LOAU  J  ) 

sC  iJ 

L  E  P  <  T  (  3  ) 

» t  W  P  ! 

CSml( i  ) 

-<5(  1  ) ; 

A 1  I  T  (  I  ) 

1 ; 

LGAn(? I 

*  C  2 : 

C  S  H  R  (  2  1 

*  » 

CSHI.  (<>> 

^^; 

CAOO(O) 

^C2: 

S  T  L  f  0  ) 

►,1)1  s 

CSHP( U  ) 

■?  a ; 

CAOO( j) 

sCd: 

A(_TT(  J) 

«-l  ; 

C  M  H  T  P  (  n 

?n  ', 

SLIT( n 

=  1  : 

I   ,  .i.KNI 


FU 
tsE 
*SK 
*SE 

i'SA 
'SA 
ASA 
<FE 
ZFi. 

*rl 

*r,H 
:*  4  H 
i!F 
tFP 
«FL 

Mil 

«pE 

*h  Y 
*SH 
<AU 
<TN 
<FL 
?  F-  E 
MF 

<  ■• ' !  I 

MN 
*FE 

.1,1 

ASH 

«.  a  n 

*SA 
'LE 
ASE 
A  rji» 

*pil 
*SE 


P  M  E  X 
T  FLA 
1^  til 
T    FLA 

vE  Ar 
VF     AC 

V  F  k  E 
tt.m  f 

TCH    F 

no   Th 

EATEP 

AT     Pfl 

ZFPO 

ASF  L 
SE.     S 

V  T  Nil 
C'JNST 

SHlF 
IFT    F 

I)   si/ 

C^fmE 
TCP  S 
TTH  S 
/FPU 
LT     iMIl 

CPF..E 
TfH    T 

V  MY 
IFT     F 
I)     PL)  I 

V  F    F  L 

avf   s 

T  f;U|l 
E  LIS 
SI  Till 
T     TMC 


T  Th 
r,   Fu 

XT  S 
d  F 1 J 
A  PO, 
A»l  . 

T  '  1 P  N 

iwsr 

TPST 

F    SH 

THA 

WFR 

.     PL 

EAO 

IHTR 

TMF 

RlirT 

TIN«j 

LAG. 

E     TU 

NT     A 

F  C  ON 
FCllN 

PET 

«Y 

MT  A 
HTRU 

6 '4  • 
I.  AT, 
HTEK 
Aii    P 
T^F 
NT     L 
S    FU 
N     I N 
prMt 


iv F  L     sTATTmFNTs. 
P     Nl)     INwFRsE. 
T ATE^FNT . 

R     INWEKSF. 


ADDMFSS. 

A  P  d  11 M  E  n  T    A  D  D  P  E  S  S  . 

A«G"M£NT    *ST'/F-*. 
ALLEST    Pp.;FH    OF    TWO 
N    UP    FQilAL    TO    SIZF. 
TS    IT. 
T11HN     ANl)    P  LAF,    EpRflP. 

HnE    FLAG  . 

ACT  PIT  nr  From  63 
PRUPFP  Pn.-R. 
Sl/F 


FLAm  HFC,  ISTFR. 
Pr,  pnlNTFP. 
i)  APr,i!MFHT  AnnRESS. 

0  APr,i|HFMT  *  N 1 1 M  h  E  R  *  . 
URN  AND  FLAG  FHPUR. 
SIZE  AN, i  fiTV  PY  64. 

P  li  P  n  I  N  T  E  P  . 

ARGUMENT  *PTTNTEP*. 
CMau'E  PF  AOORES?). 
RFGlSTfcP. 

TU  rLAli  PFGISTER. 
FuISTFR. 
IM  AnOK  FLO  HF  kCO. 

1  M  I  T  . 

P  LOOP  CONTROL. 
c.pEmfmT  rjrLn. 
i»T  T  n  U  n  f  . 
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CROTRC  3) 

i<>; 

SLIT  CO 

»o; 

CADDO) 

SC25 

SETE 

F.Un.-t, 

SETFl 

E  •  0  R  •  •  k-  i 

C  H  W  S 

64? 

Jump 

HEHR  J 

CLCC2) 

• 
• 

compc(? i 

• 

Lntr O) 

%03; 

Lt)L(  i  ) 

%\j*\ 

LOL( 3) 

*02: 

EXCHLC  3) 

*  I  C  -? ' 

End 

FFT. 

*PUSIT1HN  AODRESS  EIFLH. 
*ZERH  OUT  ADDRESS  F"!£Ln. 

*set  a  dress  FiELn  to  pointer. 

*TIJRN    UN 

*  ALL    Pt-S. 

*WlJHr>     SIZF     TO    64, 
*<iO    TO    SCRAMHLING    RiliiTlNF 

*frrhr  (no  si?E  r,TVEi. 

*SET  AGAR  ?     TO  A|  L  OnES* 
^RESTORE  APARO, 

RESTORE  AcARlt 

SFETCH  RFTDRN  ADORES*, 

* wE T ' ipix    To    CALLING    ^Por,R&H. 
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LISTED    nNI     08/  1  J/  /  1        ul  *     09  »H 


*       it 
v      x 


h  f:  r,  i  n 


X 


<;THA  ImL  1  •■  i       SECTION 
4       ******       X       X 


*        X 


X       X 

X       * 

X       X 


SIZE 

INOE 
^ODE 

lohp 

COUN 

* 

% 
* 
% 
% 
t 
% 
% 


I 

Xl 

I 
I 
T: 


EXTERNAL 
Eqi. 

Edll 
HFf'C 

Eqii 

kEi,c 
ECU 

HEf,C 

EOU 

HEGC 


M'TI  .PElfvUF.Xl 

K  S 1 1 1  ■  4  C  '■)  l 

n ; 

M  ,4 1 1 F  x  s  %  (1 0  | 

o  [  l  • »  e  ■  *  C 1  : 

i  ; 

*  L  i ) ! )  f*  s  *  r     5 

irlJU  jT  »>r  a « 
J  i 


THE    PE    MF.MnpY    L'lCATllHM    PH^nlX    CONTAINS    A    SET    <1F    STX    ImoECES 
FUh    FACH    PL    PACKED    MTU    flfit     «MR<i(SI  X    BITS    PER    INOE*.     3A    TOTAL). 
Tnr     nEF^t    NAMFp    pETr.HlNuE,,     „HI*H    FOLLOWS    mASkS    OFF  qTTS    f«Um 

pEtnhLx    ANU    STORES    ThE     uFSULl     I  >i    REGISTER    X    AND    REGISTER    fit       ThE 
PARAMETER    TO    THf    DEFINE    UFTERMIUFS    "HlCH    HITS    UF    PElNOE*    ARF    TO 
HE    KTCHEU. 


OEFINi 


FETCHlNOE X 

LOA 

SHAHR 

CLPA 

ShAhL 

LDS 

LOX 


'  STARTSahOwvIANYJ 


PEI  ii)FX: 
A  4  -  n  S 1  A  K  r ' 

ff 

s  a  ; 


*FETCH  PACKED  iNnFCEs. 

*ALIGN  TO  STARTING  MJT  POSlTTDN. 

*7ERii    REGISTER    A. 

JrEcovEr    rt  T  TS    OF    INTEREST. 

«]NnFX     Til    REG    S    ANU    DFr,     x. 


FKnM 
ANr 

STaK 
IN  A 
TlnN 

j  [n,H 

A    A 
THr 
THr 
CU  P 

I  N.- 
St q  A 

OUrU 
DEFINc 


LACM 
ft  ST 
hETOR 
TlN(, 
LL  PF 
.  RE 
LUCA 
Y  THA 
C,U  MH 
ONE  I 

using 

OATA 

MHLIM 

mLNT 


dpEkat 

AKTlNf, 
N  1  N  G     T  m 
LUC  ATT  n 
S    S  T  "1 1 1 1 
G I  S  T  E  R  s 
T  1  (1 N  S    R 
T     THFY 
L>E     Pf  G 
NPUT     Pa 
THF     nA 
SETS    LA 
G     DFAI.T 
UF     SEPT 


I  UN     i)F     ThF      ,rKAM,>Ll  Nf, 
L  DC  A  T  1  Uii     IN     '"  F      >E     Tfi 
r     PhF  VI  IUIS    C  INTF  NTS    ilF 
M.         I  N  irfHA  .j  ,L     IS     A     L)E 
TAJfcJJUSLY     JHfH     <1iIPPLTE 
X     ANU     S     C  tl  g  T  A  T    t     INDFC 

FbPLf.  1 1  'Fly.      i. if    than 

r«M  it  T  rEHATLn  hy  sin 
lb  ASbtJ  iF  0  Til  C.  inTaTn 
RAHE.TEK  RE  njIWEfl  I  S  TH 
TA  St-  T  .  THIS  DFFlNE  * 
Ri'iErf  THA  -i  1;'-  »  A  NO  IT 
Nli  .%  II  H  S  JUArtt  ^ARTTTl 
FM-FP     1^0.     uLL     ''FS     AH 


H'lilTlNE 
A  ')LSTI 
TMF_  l)F 
FTmE  fM 
I)  "JlTH 
ES  KIJH 
SFFHS  A 
PLY  PHI) 
THF  PPM 
F     NUMHF 

ILL    DNL 
CriMPLE 
ONS    AS 
E    ASSU* 


INvOLV 
NATION 
STINA1 J 
ICH  ACC 
THE  FOL 
S  T  A  R  T  I  m 
RL  PATT 
TING  Th 
PER  ENA 
H  DF  HO 
Y  HE  US 
TES  THF 
»ESCRI~ 
EU  FNA  ■ 


Es  pnii 

IN    AMD 

iin   jn 
omplis 

luring 

d  AM) 
ER'JFn 
F  TnDF 
•<LlNG 
•tS  OF 
Fi)    IN 

PtlRTI 
EU    TN 

LEn. 


T  T  n G    DATs 

TuFR    pE. 

Tuf 

HFS    THTS 

TNFORMA- 
DTSTINA- 
IN  SUCH 
CFS. 

PATTERN. 
STDRAGF 
ScpAh^l- 
(JM  OF 
THF  FFT 


OLOOPi 


I  N  T  F  R  C  H  A  m  < ,  r 

BFGTN  r^LUC^ 

STL(MUDE) 

LI  Tf LOOP) 

LIT(  INUFX  1 

RTL 

LOS 

HTL 
LOX 


♦PATTERN : 

1  »KNiU1/2»l  .  1  ; 

'  31 

t~  s  .  1  : 
•  w; 

sx»-l  ; 


THE  SYmmdL  DLOOP  MUST  Rr  LOCAL. 
*SA\/F  H(JDF  PATTERN. 
XLUUP  COUNTER, 

XPOUTTNQ  TNUEXCLOUP  COUNT). 
<iiESTTNATinN  INDFX. 
tCOMPLETE  S  TO  S  TRANSFER. 
iSlAWTlNG  INDEX. 
^COMPLETE  X  TD  X  TRAigSFEP. 
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LOR 

lda 

HTL 
LnEEl 

LljA 

SETE 

SETE1 
LOR 

STA 

HTL 

LOA 

KTL 
CROTR(mOdE) 

LOEE1 
LOA 

SETE 

SETFl 

STA 

ALITC INDEX) 

CROTLC  MODE) 

COR(MUnE) 

TXLTM(LOUP) 

RTL 
LOS 

LnR 

HTL 
LOA 
HTL 
LoEEl 

LnA 
SETE 

SETFl 

STA 

ENO 


*0(COUNT)  \ 

-4NOH1 
KiOoEi 

,:<RI 
E.UR.-U 

E.Oa.-fci 
«0(COnwt)i 

*0(cOUUT) t 

0(lMOE*>* 

SRI 

rNIMJ 

*NIJ  i(  lNnEA)  I 

«mooe» 

>bH* 

E.UP.-Ll 

E.Or.-Ej 

*O(0DU  U  )  I 

■  -1I 

*NUM(  IN0£X)| 

SPATTERS j 

.ULOOPJ 

*  S » 1 1 

«0((-nijNT  >; 
&nu  ••/?; 

*RJ 

-&NUMJ 

smUoE  i 
-  R; 
E.O-'.-ti 

F.n-.-t: 
»0  (  ("JOINT  )  S 
>  t*  : 


*FETCH  STA 
JR'JljTE. 

XSAVE  STAR 
*FAK£  END 
JENO  AROUN 

IFILL  IN  S 
HNAHLt  AL 

*FETCH  nES 
*STORE  STA 
xR0UTF  OES 
*SAVF  OEST 
*FAKE  END 
*FNn  AROUN 

% 

*FILL  IN  0 

*F.NABLE  AL 

*STQRF  OEs 
KINCRFHFNT 
^RESTORE  M 

kmOoIfy  MO 

•RETURN  TO 
XoESTINATI 
'COMPLETE 
<FETCH  STA 
XF*CHAN(iF 
*SAVE. 
tFAKF  LNO 
*ENn  AHOMN 
<FILL  In  n 
«ENARLE  AL 

iSTORE  E*XC 


RTING  DATA, 

TING  DATA, 

AROUND  ROUTE  (RIGHT). 

0  BY  IROM  pES. 

TARTING  DATA. 

L. 

TINATION  DATA. 
RTING  DATA  in  OFST» 
TINATION  DATA  rack, 
INATIUN  DATA. 
AROUND  ROUTE  (LEFT). 
0  BY  UNijM  pEs» 

ESTINATION  -)ATA. 
L. 

T.  OATA  IN  (;TART  LO'". 

index  (■-Luop  count), 
ooe  pattern  (almost). 

Dt  PATTERN, 

LOCAL  DLOnPflNC  LOOP> 
ON  INDFX. 
S  TO  S  TRANSFER. 
RT  AND  OEST.  DATA. 
START  AND  DEST.  DATA. 

AROUND  ROUTF  Or  4NUM/? 
D  BY  4NUM  PES. 
ESTINATION  DATA. 
L. 

hAnrEo   start    And   dEsT. 


FOR  THF  SMALLER  DAT/.  6FTS  (*>■'*  Akd)  S-iALLFR)  'SCRAMBLING  CAN  HE 
ACtUmRLT  ShE{)  vITm  UNly  HNE  ACCFSS  TO  ONT  K<\m    Of  sTOpA(;F,   Tu|F 
INTERCHANGE  OPEpAUHN  CAN  BE   UflNE  USING  ONLY  REGISTER  TpANSrFRS 
AN-,  RUljTINvi,   E/CHAmuE  IS  A  OtFTME  «hICM  PERFORMS  ONE  TnTFRCHANgF 
ON  UATA  FOUND  In  Kt,ISTFK  A,   THE  DISTANCE  SEPARATING  ThF  TWO  wORD.S 
TU  BE  InTERCHANgFD  IS  T  nt  UNLY  PARAMETER.   A  proper  mudf  pattern 
MUcT  BE  SUPPLIED  IN  THE  CD   nlUE  REGISTER. 


D E F  T  N r   EXCHANGE  t  Nil  i  J  » 

LnR  <*A; 

RTL  -,uu  <; 

LDEEt  ».modf; 

lds  «.a; 

Li>A  VR? 

RTL  AS. ->c  1 1  DM} 

CSHL(MI)OE)  C.NIJMl 

LDEF1  htU'jUti 

Lda  »ii i  it m 


<<nATA  TO  REGISTER  R. 

M<OUTF  STARTING  DATA, 

*SELFCT  ELEMENTS  TU  RF  EXCHANGE' 

^PERFORM  THE 

^EXCHANGE . 

«ROl|TE    RACK    TO    ThF    STArT, 

^POSITION  MUDF  PATTERN  TO 

* F  ILL  IN  DESTINATION  DATA, 


% 

HERE 


t  E  N  T  R  Y  1  :  « 

CTSHT(ST/F)  <SO»X«; 

CTSRT(ST/Fi  VV#A)6S 

CTSRT(SIZE)  bd.xi^; 

CTSRKSI/F)  Sf.xisn; 
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Xfil 

M6« 

X32» 

X6a« 

^128 
X256 
X512 
X10? 

X2oa 
i 
x 
% 


\ 
i 
i 

4  I 

8t 


CTSHT(SIZE) 

CTSHT(ST^F) 

CTSrT(SIzE) 

CTSRT(SIZE> 

CTSRT<SIZE) 

CTSBT(SUE) 

CLC(2) 

CC1MPC(2) 

LDL(O) 

LDL(1 ) 

LnL( 3) 

LXCHLO) 

Jump 
Jump 
j  u  m  p 
jump 
Jump 
jump 
Jump 
jump 
Jump 


<>6.  <  1  2  C3  » 

53*X1024) 

52#X?04H> 
Si  #440*61 

; 

j 

SD  J  5 
SO 4  I 

M)2, 

T  icr? 

a«  J 

A  1  6  » 

ft  32* 

u64  : 
A128' 

a2b(S  5 
Abl?  » 
a1 o?a J 
A204HI 


tERRHRt  DATA  SlZT  Tun  LA*r,E. 
*SET  ACAR?  Tl)  ALL  ONES. 
*RESTDkE  ACARO. 
^RESTORE  ACArtl. 
XFLTCM  RETURN  AonRESs. 

KRETHRN* 


this  portion  uf  th^  pkogkah  reorders  <scramrles>  The  data  in 
surm  a  way  that  the  results  uf  the  toansfurm  routine  will  come  out 
st°ufntially  stored  in  or  mfmury.  df'ly  one  of  the  follo-mng  ten 
strhuns  uf  codf  "ill  he  executed  depending  on  the  sl*e  of  the  data 

SEtS.   THE  SCRAMHLlNG  PROCESS  IS  PURELY  A  DATA  TRANSFER  PoOrlEM 

AN<x  TS  ACCOMPLISHED   HJICKLWA  FRACTION  OF  THE  TRANSFORM  TlMF).   THF. 

pRTMARY  OptRATInNS  ARE  MLHUrtv  ArcESS  A,Nn  HOuTlNr,,   T  hE  oThEr  IN- 

S'dUCTION  SET  UP  INDEXING  Ami)  L"DP  CONTROL. 


i409MI    FFTCHlNOE*  5f»*S 

LIT(LUOP)  1»J1.1» 

LIT(IndEx)  *b*l 

LUOPIj     MTL  *SM! 

LDS  *RS 

RTL  M.-l j 

Lox  &R! 

LOR  *0( COUNT) 

RTL  O(LUOP); 

LDA  *0(Cr»UNT) 

STR  *0( cDUNT ) 

RTL  *A»0(I*  )F. 

STR  *0(rOUNT^ 

ALIT( IN  >FX)  s-1  ! 

TXLTM(LODP>  »ludp1; 

R  T  I.  i  S  .  1  J 

Los  « R ; 

LOR  i 0 ( n  O I » N  T  ) 

RTL  12: 

STR  *0< COUNT  > 

ALIT(COIINT)  r*3j 

TXLTMfC'IIMT  )  ,A409h; 

JUMP  FFTl! 

A2046  t  t  FETCHINDEX  S 1  J  6  : 

LIT(MIJUF)  .100  0  0000 

lNTERCHAMr,r  32; 

ALITCCOUNT)  =31! 

TXLTM(COUriT)  »A2f)«a? 


J 


<0E 

<L0 

srU 

•<dE 
*CU 
*ST 

*c'J 

*FL 
•<H0 
*FE 

1ST 
*R0 
''  *ST 

XIN 
*rE 

adE 

*co 

J     *fE 

*R0 

I  *ST 

%\H 
*r,o 
*JU 
SUE 

ooopooooo 

xoE 

*IN 


ST 

l)E 

oE 

>EN 


FINE  ( 

np   cnu 
ut  I  nc, 

STINAT 
MPLETE 
ARTING 
MPl  ETE 
TCH  ST 
UTF« 
TCH  DE 
ORF 

UTE 
ORE 
CRF' 
TURN  T 
STINAT 
MPLtTE 
TCH  ST 
DTE  EN 
ORF  EX 

CRF.'EN 

RACK 
MP  TD 
FINE  ( 
OOOOOl 
FINE  ( 
CRFMEN 
HACK 


FETCH 

NTER. 

INDEX 

ION  I 

S  TO 

INDF 

X  TO 

ARTIN 

END  A 

STINA 

ARTIN 

STINA 

ST.  n 

T  IND 

0  LDO 

ION  I 

S  TO 

ART  A 

D  ARO 
CHANG 

T  COU 
FOR  N 
TRANS 
FETCH 

HI 

LARGE 
T  COU 
FOR  N 


FROM  PEINOFX). 


(-LU 
NDFX. 

S  TR 
X. 

X  TR 
G  OAT 
ROUND 
TION 
G  OAT 
TIOM 
ATA  I 
EX  (■ 
P.lf  IN 
NDEX. 

S  TR 
NO  OE 
UNO  b 
ED  ST 
NT. 
EXT  D 
FORM 

FROM 


Dp  fUUNT), 

ansfEr. 
ansfEr. 

A. 
RY  64 

OATA. 

A  In  oFST. 
data  rack. 
n  start  loc. 
-loop  cuunt), 
c  loup  count) 

AnsfEr. 

ST.  DATA. 

ART  ANO  DEST. 

ATA  SFT. 
ROUTlNEt 
PEINOF*). 


). 

NT. 

EXT  DATA  SFT. 
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A  1  0  2  4  I  I 


L00P2I 


A512I I 


L00P3t 


A256| i 


L00pat 


Al2«l  I 


loopsj 


jump 

FETCHlNOEX 

LIT(MODE) 

iNTERCHANr.F 

LIT(LUOP) 

LDA 

LIT(MUDE) 

t-xCHANGE 

SETE 

SETEl 

STA 

ALIT(CI)IJNT) 

TXlTMCLOnP} 

TXLT(COUNT) 

JUMP 

FETCMINUFX 

LIT(MUI)E) 

iNTFRCHAN^r 

LIT(LOOP) 

LDA 

LIT(MOl)E) 

EXCHANGE 

SETE 

SFTF1 

STA 

ALIT(COUNT) 
TXLTMCLHUP> 
TXLT(COUNT) 
JUMP 

FETCHlNOEX 
LITCMODE) 

INTERCHANGF 

LIT(LUnP) 

LOA 

LIT(MOOF) 

EXTHAN'iF 

LIT(MQOE) 

EXCHANGE 

SFTE 

SETEl 

STA 

ALIT(COUNT) 

TXLTMCUHJPi 

TXLT  (  CHIINT  ^ 

JUMP 

FETCHINOEX 

RTL 

LOS 

LOA 

LIT(MUi»F  ) 

EXCHANGE 

SETE 

SETEl 
STA 

LIT(LOUP) 
LOA 

LIT(MUDE) 

E  X  C  H  A  N  r,  r 
LI  T(  -iUDD 
EXCHANGE 

SFTE 


FFTU  XJUMP  TO  TRANSFORM  ROUTINE. 

4t>»*>  XDEFINE  (FETCH  FROM  PElNnFX). 

■1000004000020 OOOl OOOOO'M 

16|  ^DEFINE  (LARr.E). 

1  » 1 5 • 0  *         *LOOP  COUNTER. 

O(CnUNT)!        *FETCH  DATA. 

■  iU/fftoOOUOO'B' 


toEriNt.  . 

«ENARLL  ALL* 
* 

*STORE  RESULT. 

^INCREMENT  COUNT. 

XGU  RACK  FOR  NEXT  ROW. 

*GO  RACK  FOR  NEXT  OATA  SFT. 

*JUHP  TO  TRANSFORM  ROUTINE.' 

UnEFINE  (FETCH  FROM  PETNnFX). 
*  10  02  00a 01  00200'! 01  0  0200  IB  I 
rti  *DEFINE  (LARGE). 

1 »f #01  *LUOP  COUNTER. 

OlCOUNMJ        XTETCH  OATA. 


E.UR."El 
E.U'v.-t  * 
(KCOUNT  )  l 
*\i 

■L00P2I 
XCUUNT.  A10?/4  : 
FFT1  i 
39IOJ 


24  J 

F.Uw.-LJ 
F.U^.-El 
OiCUiNl  )  ; 

si; 

.LGllPJJ 

$CUUNT» A^l?l 
FFTl  i 

33IM 


*0EFI,ME.    " 

*ENAR|.t    ALL. 

% 

iSTORE    RESULT. 

MNCWFMENT    COUNT. 

*G0    RACK    Ff)R    NEXT    RllW. 

*G0  PACK  FOR  NEXT  DATA  SFT. 

<JUMP     TO    TRANSFUWM    RriilTlNE. 

'fcuEFINE    (FETCH    F»nM    pElNOF.X). 
■104  21 04^104  21 0«2lO«2lOl«l 
a;  XOEFINE    (LARGE). 

t.J.Oi  *LUOP    COUNTER. 

0CCHUN1  )  I  *FETCH    OATA. 

*  3  h  o  (  a  1  7 ,  i  J »,  o  }  8  ; 
2«J  KfiEFTNt. 

a  f  f  ^  00  00  Ot>  ot  77  a  On  » 8* 


E.U7."LJ 

F.i I  .'.-LI 

o( coun r  i ; 
=  i; 

•  LUfiPa ; 

ftCUiluT.  A/»b6  1 
F  K  T 1  : 
?f'.t  : 

*  S .  II 

*h; 

*o(  rmiNT ) ; 


SOEFTNE, 

<F.NAhll    ALL. 

t 

*  ST  ORE    RESULT. 

%\HC.»r  ME  NT    CO'iNT, 

%r,[)    BACK    FUR    nFaT    ROW. 

*G0    HACK    FUR    NEXT    UATA    SFT. 

*JUMP    Til    TRANSFORM    »0UTlNF.. 

<i>EFIwE    (FETCH    FROM    pFTNnFX). 

*KOUTF     INDEX. 


^COMPLLTF    S    Til 
*FETCH    OATA. 


S     TRANSFER. 


=  5  2S2S2«S2t>?,S?5?,3252e>2t)l«l 

1  ;  '.OLFHL. 

F.Hy.-LI  *ENARLt    ALL, 

E  .  U  *  .  •  1 1  * 

tO(cOUNT)!  *ST0RF    RESULT. 

1  . 1  »o;  <t  UOP    CflUNTER. 

0(CHUNF  )  I        *FETCH  OATA. 

s3la(SJl4o3iaiBI 

30:  inEFINE. 

s?41ff41^n01703770360*Bj 

1?;  SOF.FTNK. 

L.U'.-t.s  iFiMARLL    ALL. 


-35- 


AMI  i 


A32»  » 


A16>> 


ABl  I 


SETE1 
STA 

ALIT(COUNT) 

TXLTM(LOOP) 

TxLT(CtHJNT) 

JUMP 

LDA 

LlT(MOnE) 

Exchange 

LIT(MOOE) 
EXCHANGE 

LIT(MUOE) 

Exchange 

SETE 

SETE1 

STA 

TXLTM(COUNT) 

JUMP 

LDA 

LIT(MODE) 

Exchange 

LlT(MOnE) 

EXCHANGE 

SETE 

SETE1 

STA 

txltm(count) 
jump 

LDA 

LIT(MODE) 

EXCHANGE 

LIT(MUDE) 

EXCHANGE 

SETE 

SETE1 

STA 

TXLTM(COUNT) 

JUMP 

LDA 

LIT(MODE) 

EXCHANGE 

SETE 

SETE1 

STA 

TXLTM(COUNT) 


E.OR.- 

0(C°UN 

■  II 

.L0DP5 

SCUiJNT 

FFT1  I 

0(CfJUN 

•25252 

311 

■63146 

M» 

■36000 

4j 

E.OR." 

E.Op.- 

occoun 

.  A64j 

FFUl 

0(CfiUN 
•52525 
15j 

■31400 

6» 

E.UR.- 

e.o*.- 

0(C0UN 

.A32> 

FFT1* 

0(CUUN 

■25200 

7} 

■60140 

2> 

E.UP.- 

E 1 0  s  • " 

0(cnnN 

,A16J 
TFTl  » 
0(C  1UN 
■50120 

3; 

E.Utf." 
E .  0  n .  - 

0(cf'dN 
»A*t 


I 

»AUH 

T)  I 
52525 

ooooo 

17000 

L» 
LI 

1  H 


T)l 
OOOOO 

U600 

LI 
El 
T)l 


T)J 

12500 

30060 

El 
Li 

T)l 


24050 

LI 
El 

Til 


*STOrE  RE 

XINCRCMEN 
XftO  BACK 

I    XGO  back 

*JUMP  TO 

XFETCH  DA 
?I8I 

XOEFINE, 
1463141*1 

XqEFINE. 
0740003601b' 

XDEFlNE, 

XENABLE  A 

X 

*STORE  RE 
XGO  BACK 
XJUMP  TO 
XFETCH  DA 
12525?I8I 

*o£finl, 

0630003141b! 

xdefinl. 
xenable  a 

% 

«TORE  RE 
Xf,0  RACK 
XJUMP  TO 
XFETCH  DA 

052400252ltt| 
XOEFINE. 

1403006014«8I 
XnEFlNE. 

XENABLE  A 
X 

XSTQRF  RE 
Xf,U  BACK 
XjUMP  TO 
XFETCH  DA 
1202405012181 
XOEFINE. 
XENABLE  A 
% 

USTORE  RE 
XGO  BACK 


SuLT. 

T    COUNT. 

FOR    OTHER    RflW. 

FOR  NEXT  DATA  |FT. 

TRANSFORM  RQUTlNE* 

TA. 


LL. 

SULT. 

FOH  NEXT  ()ATA  SET. 

TRANSFORM  ROUTINE* 
TA. 


LL. 

SULT. 

FOR  NEXT  2  DATA  SETS. 

TRANSFORM  ROUTINE. 

TA. 


LL. 

SULT. 

FOR  NEXT  4  DATA  SETS. 

TRANSFORM  ROUTINE. 

TA. 


LL. 

SULT. 

FOR  NEXT  »  DATA  SETS 


END 


HERE. 


-36- 


ilSTEO  nfy{  06/13/71       Ait  09U 


%      i 
%   % 

%      i 

sm« 

MGDEs 
LCOPj 

COUNT  t 


BEGIN 

*  ?.  * 


x 


x  *  x 

TRANSFORM 


%       %       %       %       % 
EXTERNAL 

KE6C 

EQU 

Krr.c 

R?GC 

EQi.J 

REGC 

THE  FOLLOWING 


S   8 

SEC 


ION 


%       % 
S I W  | 

01 

SSIZE**C0* 
1  J 

SMlJf)E«*Cl  * 
?'  I 

RLOOP«*C' 
3? 

*cnijNTBsrii 
OEFImE    is    a 


X   % 


%  % 


%     '& 


STANDARD  transfer  operation  used  tn 

Trtr  fU   irIEn  TRANSFORM  ROUTINE.   THE  dARTTAl  RESULTS  OF  EAcH  ITERA- 
TirN  ARE  USED  In  THE  PL  .mHLRE  THEY  *ERE  CUMPIlTED  AS  WELL  AS  IN  A 
PL  SRME  OlSTANCF  AWAY  FnR  1 hL  NF'T  ITERATION.   THE  FUNCTION  nF  THIS 
ntp-lNt  IS  TO  MOwE  THESE  PARTIAL  RE SULT S ( FOUND  IN  REGISTER  A)  TO  THF 
pHopER  pES  uoIST  IS  A  NumhEP  TELLlNf,  HOvg  FAR).   T  hF  pROpFR  hOqE 

pATTFRN  MUST  BE  SUPPLIED  IN  THE  CU  REGI STER  .MODE ( AC  A R 1  )  .   NnU  ThAT 
RL-IsTER  MODE  (ACAR1)  IS  NUT  CHANGED  BY  THIS  DEFINE, 


DEFINr 


TRANSFER 

RTL 

COMPC(MODE) 

LDEE1 

LDS 

RTL 

COMPC(MW>E> 

L0EE1 

LDS 
LDA 
SETE 
SETE1 


*DIST»» 
S  A  ,  X  D I  S  T  I 

I 
$ mUoEi 
sr; 

SA»-4DIST> 

t 
RMUnE} 

<a; 
$rj 

E.UR.-El 
E.Ow,-Ei*»i 


*ROUTF  REGISTER  A  RlfiMT. 

^COMPLIMENT  MODE  PATTERN. 

IENARLE  PES. 

<LOAD  REGISTER  S  IN  nESlREO  PES. 

*ROUTE  RFGISTER  A  LEFT. 

^COMPLIMENT  MODE  PATTFPN. 

XENARI.E  PES. 

XLOAD  REGISTER  S  IN  oFSlRED  PES. 

*LOAD  REGISTER  A  IN  SAME  PES. 

^ENABLE  ALL  PES. 


PAGE 


I 


THE  following  define  pehfoms  THE  COMPLEv  multiplication 

OPrRATlON  ON  Two  UPERANdS.   EACH  OPERAND  IS  ON  wORO  COMpQsEn  Op 
T*n  32-BIT  FLOATING  PUlNT  NUMBERS*  THE  OUTER  OF  wHlCH  Is  THE  REAL 
pApT  AND  THE  INNER  IS  THE  IMAGINARY  PART.   THE  TwO  OPERANDS  MUST 
RE  In  THE  A  AND  S  REGISTERS  IN  ALL  PESt   THE  RESULT  IS  |_EfT  IN 
RE-.IsTER  A.   THE  RESULTS  IN  TErvS  OR  RKREAL  PART  OF  FIRST  OPERAND) 
RiJ.Il.ANn  12  ARF  AS  FOLLOWS, 

R«R1*R2-I  1  *\'i 
l«Rl*I2+R2*Il 


DEFINr 


COMPLEXMULT 

MLRN 

LOR 

SETE 

CHSA 

SETE 
LOS 
LOS 
LDA 

SWAPA 
M(_RN 

LDB 


it  W  H  "l  I  * 
%Sl 

*  WHQ  I 

E.Aij  ».-E 
l 

E.O..-L: 

SA; 

SH{ 

I 

SRS 
SSI 


^MULTIPLY  REALS  AND  IMAGINARIES. 
*WEL0AU  SFCOND  OPERAND. 

*CHANr,E  SIGN  OF  INNEp  (I**2*-l). 

* 

triNE  OPERAND  TO  R  REGISTER. 

*SAVE  PARTIAL  RESuLT. 

*FIRST  OPERAND  Tn  A  REGISTER. 

XIMAG  TO  OUTER*  REAL  To  INNER. 

XMULTlPLY  R1*T2  AND  R2*Il. 

^POSITION  PARTIAL  RESULTS 


-37- 


IN  HE 
ATr  m 
SHnUL 
G I  M  E  N 

MUi  Tl 
CUmST 

TMr  p 

TlrN. 

ALl  UW 
TlnN 

ON  CO 


SWAPX 
AON 
THE  FU 

gIstEr 

Emory 
o  re  e 

ITERA 
THE  BA 
RLlCAT 
ANTS  ( 
REVIOU 

THE 
S  HALF 
TO  THE 
MPLEX 


OEFINr  M 
L 
C 
L 
L 
A 
S 
I 
S 

s 

T 
ANn  CO 
CUf'STA 
wAc  TA 
OF  AnG 
OF  ThO 
FIpST 

1 
THFM  I 
OEc InE 
STnKEO 
TU  BE 
RErUfH 


ULTAD 

DA 

OMPLE 

OR 

OS 

ON 

TA 

OA 

RN 

TA 

HE  CU 

sines 

NTS  R 
KEN  0 
LES  I 
SE  OF 
OlJADR 
HE  FU 
NTO  T 
ARE 
CONS 
THE  A 
ED    IN 


LLO*I 
S  An 
LOCAT 
QUAl 
TION. 

sic  n 

ION  r 
W  S). 
S  ITr 

SYMET 
THE 
OTHF 

MIJLTI 

OSUR 
XMULT 


I 
SB) 
N<i  DEFI 
D  PRODU 
IUNS. 
TU  HALF 

THE  S 
PERATIO 
ETwEEN 

THIS 
RATION 
RY  OF  T 
OPEHAT1 
R  HALF. 
PLY  UNE 


##> 

NE  HEUUI 
CES  T  ,U 
THE  OFFS 
THE  NUM 
IGNIFICA 
N  OF  THE 

the  resij 
new   «uAn 

TU    PRUOU 

HE    CALCtj 

ONS    Tl)    R 

THUS    A 

AIJU     AND 


XFOR  FINAL  ADDITION, 


»ES  PR 
RESULT 
FT  USE 
PER  OF 
NCE  OF 

FAST 
LTS  Of 
TITY  I 
CE  RES 
LATION 
F  ELIM 
STAND 
ONE  S 


OPER  CON 

S  AND  $T 

0  Tn  FIN 

ROWS  BE 

THIS  OE 

FOURIER 

ONE  ITE 

s  Then  a 

ULTS  FOR 
S  AS  «FL 

INATEO  R 
ARD  OEFI 
URTRACT 


STANTI  TO  R 
OREs  ThCm  I 
0  TMC  FIRIT 
ING  ORCRATf 
FINE  FOLLOW 
TRANSFORM  I 
RATION  AND 
ODEO  TO  THf 
THE  PRESEn 
L  AS  THAT  0 
Y  ADOINQ  ON 
NE  WHICH  AC 
AND  PRODUCE 


t   Fnu 

N  ApR 
ORCR 

0  On 

IN  B 

1  A  C 
A  SET 

rEsu 

T  ITE 
F  The 
E  SlJB 
CnNPL 
S  Two 


nd 
OprI- 

ANO 
IN  A 
RIEF. 
OMPLE* 

OF 
LTS  OF 
RA- 

w  S 
TRAC- 
ISHES 

RES* 


NSTAn 

RATh 

EoUlP 

F  THE 

N  ANY 

ANGL 

ANT  S 

LLOWT 

HE  PR 

LEFT 

TANTs 

ODREs 

THE 


* 
& 

0 
*, 
S 
0 

s 
$ 

& 

(« 
TH 
TO 
RE 
T 

ES  IN 
I  WES 
NG  OE 
OPER 
IN  RE 

ARE 
S  UF 
CURHF 


TS 
EH 
EU 
IK 
UF 


OFFSE 
OFFSE 
UFFSE 
(CniJN 

A; 

R» 
(CHUN 

R; 
SJ 

OFFSE 
S)  F 

AN  EX 

BE  S 

GULAR 

HE  FO 

THE 
AND  C 
FINE 
SECON 
GISTE 
NOT  U 
THF  M 
NT  TT 


Tla 

I  CCUil 
KCUU 

I )) 


r>» 


T(CUU 
UR  TH 
HONLig 
TOHEf) 
iTv  A 
UR  'JU 
FIRST 
USINE 
LUAUS 
U  >UA 
h  S  A 
ISTUR 
EMOKY 
ERATI 


NT)|«fETch  first  opErano, 

NT)»*MULTTPLY  BY  CONSTANT 

XFETCH  SECOND  QPERANn 

*SAVE  PARTIAL  RESULT, 

*ADD  ?ND  OPERANO  TO  P 

XSTORE  ONE  RESULT, 

X?No  OPERAND  TO  RfGls 

tSUBTRACT  PARTIAL  RES 

NT)»tfl   ISTORE  THE  OTHER 

TS  PROGRAM  ARE  REpRESENTe 

TIALS.   TO  DECREASE  THE  n 

AS  PART  OF  THE  PROGRAM  a 

MO  SYMETRY*  THE  sines  or 

AORANTS  ARE  SIMPLE  TRANSF 
OUADRANT,   FOR  THIS  REAs 
S  ARE  STORED. 

A  SET  OF  CONSTANTS  ANO  T 
nRANT  VALUES.  ThE  RESULT 
ND  THE  MEMURY  LOCATIONS  0 
RED.   THE  INPUT  PARAMETER 

LOCATION  CONTAINING  THE 
ON.   THE  TRANSFORMATION  j 


(W)  IN  RQS. 


ARTIAL  RES. 

TER  A. 

ULT. 

RESULT, 
n  As  SINES 

UMBER  OF 
DvANTAGE 
COSINES 

ormation 

ON  ONLY 

RANSFORMS 
S  Or  ThIS 
F  THF 

TS  ASSUMED 
CONSTANTS 
S  DEFINED! 


KEAL(UUAO    ?)»-lMAGINANY(UUAO     1) 
IMAGINARY(QIIA0    2  )  »PE  AL(  vJUAn     1) 


DEFINr       CHANGEUtlAI) 
LDA 
SWAPA 
SETE1 
CHSA 
SETE1 
LOS 


<,CUMHLLXNUM1 
ICCUMPLEXNUMI 

; 

E#ANO."EI 
E.UR.-Ej 


XFETCH  CONSTANTS. 

KSWAP  HEAL  AND  IMAGINARY  PARTS* 

% 

*CHAMGE  SIGN  OF  REAL  PART. 

I 

*LEAVE  RESULT  IN  REGISTER  S. 


THE  FULLOWTNb  DEFINE  RESLTS  THE  TRUE-FALSE  FLIP-FLOP  (TFFF) 
WHtCH  IS  USED  TO  COWTRUL  A  SMALL  LOOP. 


OEFINr   RESETFF  = 

CACRB  0  **l  XRESET  ACR  HEfllSTEH  -R I T  0  CTFFF) 

I 

*  THE  FOLLOWING  DEFINE  SETS  ThE  TpUE-fAlSE  FLIP"FLOP  wrICH  15 

*  USrU  TU  CONTROL  A  SMALL  LOUo. 

OEFINr   SETFF  a 


-38- 


CA 
AT 
MUsT  BE 
OF  ElGH 
SUrCEED 
IN  EACH 
Sk-T  IS 

cunE  is 

THaT  SE 
CUmPLET 
FUcM  IS 
EArH 

qata 


ST 
SE 


X 

% 
% 

% 
% 
% 
% 
I 
t 
% 
% 
% 
% 

EI 
FFTHENtRY]  J 
SE 
SE 
CH 
CA 
LD 
SL 
CA 

F8»        LO 
LI 

HT 
CO 
Lfl 

Ln 

HT 
CH 
CO 
LO 
LD 
LO 

SE 

SE 
An 
LI 
RT 
CO 

Ln 

LO 
HT 
CH 
CO 
LO 
LO 
Ln 
LI 

Ln 
sw 
SE 
CH 
SE 
SE 
An 
LI 
TR 
ML 
LI 


CR8 

THIS 
GIN. 
T  DAT 
ING  S 

ITER 
USED 

LABE 
CTION 
ED  WH 

nUNE 
EP  IS 
T(S> 


POIM 

THF 

A  Pni 

TEP  T 

ATInN 

FOR  T 

LED  B 

(ONE 

ICH  C 

AMn 

CHMP 

RElNG 


LL 
i 

TE 

TE1 

WS 

CRB 

L(SIZE) 

IT(COIINT) 

DOCCOUNT  ) 

A 
T(MODE) 

L 

MPCCMRDE) 

EE1 

S 

L 

SA 

MPC(MODE^ 

EE1 

S 

A 

TE 
TE1 
N 
T(MOOE) 

L 

MPC(MnOE) 

EE1 

S 

L 

SA 

MPC( MnOEl 

EE1 

S 

A 

T(HODE) 

EE1 

APA 

TE1 

SA 

TE 

TE1 

N 

TCMOHE) 

ANSFER 

RN 

T(MOOE) 


200 
T  SCRA 
rRANSF 
NTS  CO 
WICE  T 
.  THI 
HE  rIN 
Y  THE 

STEP 
UHRFSP 
THE  RE 
LETELY 

TRANS 

16» 


IB  ##l 

MBLING 
ORM  AL 
NTAlNE 
HE  PRE 
S  PHOC 
AL  STE 
NUMBER 
OF  THE 
ONDS  T 
MAIN  C 
INDEP 
FORMFil 


HAS 
GORIT 
D  IN 

vious 

ESS  C 

•»•   T 

OF  D 

TRAN 
U  THE 
UDE  I 
ENT  0 


XSET 
BEEN  C 
HM  BEG 
EIGHT 

NUHBE 
ONTlNU 
HE  BEG 
ATA  PO 
SFURM) 

SIZE 
S  SKIP 
F  THE 


ACR  R 
nMPLE 
INS  R 
CONSE 
R  OF 
FS  UN 
INING 
TNTS 
,   WW 
OF  TH 
PED. 
SIZE 


EGISTER 
TED  AND 
Y  WORKI 
CUTIVE 
DATA  PO 
TIL  THE 

OF  EAC 
IN  A  GR 
EN  THE 
E  DATA 

IT  IS 
OR  THE 


THE 
NG  WT 
PES. 
INTs 

ENTI 
H  SEC 
OUP  U 
STEP 
SET  T 
IMPpR 
NUMq£ 


0  (TFFF). 
REAL  WORK 
TH  GROUPS 

IN  EACH 
ARE  USED 
RE  DATA 
TlON  OF 
SEO  IN 
HAS  BEEN 
u|E  TRANS- 

tant  that 

R  OF  THE 


*align  to  b  pe  bounory. 


E.Ow.-Lj 
E.Oq.-Ei 
325 

211  IBI 
SOl  j 
■  01 

ssize; 

O(CQUNT) 
.1252525 

*A»  1  J 
J 

RHODE* 
SRI 

SA.-l  | 

J 
I 

SMUnE»* 
*  A ; 
SR) 
E»UR»"Ei 

E  .Ow.-E» 

SSI 

»14631<*6 

SA,?; 
I 

SMODEJ 

SRI 

SA ,-2  ; 


irtODEI 
*.  a  ; 

jH  ! 
.5252525 

SMOOEJ 
I 

E.am  >.-r 

i 
E.JR."Ei 
E.Q4.-E: 
«s; 

*lf03f>07 

SINI 
.5252525 


*ENARLE  ALL  PES. 

X 

15 

i 

% 

% 


25252525 

% 
% 
% 
% 
% 
% 
% 
% 
t 

« 

% 

i, 

31463146 
% 
% 

% 
% 

% 
% 

% 
% 
* 
% 

2*252525 
i 
r. 

I  % 

* 

% 

% 

4l7n3607 
% 
% 

25252525 


CHANG 

TURN 

FETCH 

RESET 

SET  T 

LOAD 

2525? 

ROUTE 
CUMPL 

enabl 

LOAn 

ROUTE 

CHANG 

COMPL 

ENA3L 

LOAn 

LOAn 

ENARL 

PARTI 
31463 

ROUTE 
CUMPI. 

FNARL 

LOAn 

ROUTE 

CHANG 

COMPL 

ENABL 

LOAD 

LOAn 

25252 


E  TO 

OFF  U 

SIZE 

COJJN 

U  STA 

data. 

52lB» 

REGI 
IMENT 
E  PES 
REGIS 

REGl 
E  SIG 
IMENT 
E  PES 
REGIS 
REGIS 

E  ALL 


32-BIT  HODE. 
NDERFLOW. 

INFO. 
T. 
RTING  POINTER. 


STER  A  RIGHT* 
MODE  PATTERN. 

TER  S  IN  DESIRED  PES. 

STER  A  LEfT. 

N  OF  A  IN  DFSTRED  PES 

MODE  PATTERN. 
* 

TER  S  IN  QEslRED  PES. 
TER  A  IN  SAmE  pts. 

PES. 


AL  RESULT. 
14|fl| 

REGISTER  A 
IMENT  MODE 
E  PES. 
REGISTER  S 

REGISTER  A 
E  SIGN  OF  A 
IMENT  MODE 
E  PES. 
REGISTER  S 
KEiISTER  a 
5H| 


RIGHT. 
PATTERN. 

IN  DESIRED  PES. 

LEfT. 

IN  DESIRED  PES 
PATTERN. 

IN  DESIRED  PES. 
IN  SAmE  pEs. 


FAKE 
SINES 

A  mo 

COSINES. 

PARTIAL  RESULT. 
/ll70360l«l 
ROUTING  DEFINE. 
MULTIPLY  BY  CONSTANT. 
25?525«8» 


•39- 


NE8l 
Fl6l 


NE16I 
F32I 


NE32: 


F64I 


LDEE1 
LDR 
SWAPA 
SETE1 

CHSA 

SETEl 

ADM 

LIT(M 

LnEEl 

SWAPA 

cshr( 
loei 

SETE 
CHSA 

SETE 

SETEl 

AON 

STA 

TXLTM 

CTSRF 

LIT(2 

JUMP 

SLIT( 

CAOD( 

LIT(M 

LOA 

TRANS 

STS 

Lns 

LnA 

COMPL 

AON 

STA 

TXLTM 

CTSBF 

LIT(2 

JUMP 

SLIT( 

CADO( 

LIT(M 

LOA 

TRANS 

STS 

LOS 

LDA 

COMPL 

AON 

STA 

TXLTM 

CTSBF 
LIT(2 
JUHP 
SLIT( 

cado( 

LITCM 

LDA 

RTL 

LOS 

LDEEl 

LQS 

LOA 


*MOUE! 

$  A  I 


F.ANO.-Ei 


l 


UOE) 


MORE) 


(COUNT) 

(SIZE) 

) 

COUNT) 
COUNT) 
ODE) 

FER 


EXMULT 


(COUNT) 

(SIZE) 

) 

COUNT) 
COUNT) 

ODE) 
FE» 


EXMULT 


(COUNT) 

(SIZE) 

) 

COUNT) 
COUNT) 
ODF) 


E.OR. 
SR> 

•  2104 

smUoe 

J 

1* 

*Mqoe 

E.ANO 

1 
E.OR. 
E.OR. 
SSI 
0(COU 

•  F«» 
36.  NE 
PAIR( 
OUT) 

■  01 

8SI7E 

»1774 

o(cnu 

0(COU 
SA» 
SIN+I 
SIN  +  1 

0(COU 
0(COU 

•  F16> 
35.  mE 

PAIR( 
OUT  I 

tSIzE 

•  1777 
0  (  C  0  U 
lb? 

o(cnu 
sa  ; 
SIN  +  2 
SIN*2 
0(COU 
0(COU 
.F32J 
34,  mE 

PAIR( 
OUT* 

ssizt 

*\777 

0(COU 
SA.32 

sr; 

smuqe 

&aj 

*R  J 


21042 

I 


I 

.-E» 

-Li 

-Et 


NT)| 

s; 

4.0.4 


I 

00776 
NT  )  I 

NT)I 

I 

I 

NT)| 

NT)I 

16, 

8  .  o  »  8 


I 

77000 
Nl  U 

NT  )  J 

1 
I 

Nl  )  I 
NT)  I 

32, 

16.0. 


I 

7(!7( 

NT)I 

I 

» 


« 

X 

X 

I 
I 
I 
I 

104210 
I 
I 
* 
I 
% 
X 
% 
% 
X 

X 
t 
X 

.0)1  X 

% 

% 

X 

003770 
% 

X 

X 

X 

% 
% 
X 
% 
i 
% 

.0))  X 
X 
X 
X 

003777 
X 
X 

I 

X 

* 
I 

X 

X 
X 
X 

16.0)1 
X 
X 
X 

77«000 
X 
X 
X 
X 

% 

X 


FAKE 

CUMPLCX 
MULTIpLy 

BY 

TAKING 
ADVANTAGE 
421042lAf 
OF 
THE 
SfMETRY 

or 

SINES 

ANO 

COSINES, 


PARTIAL 
ENO  OF 
ITERATE 
GO  ON  I 
SET  UP 
QUIT  IF 
ZERO  CO 
START  W 
0177*00 
RESULTS 
ROUTING 
SAVE  UT 
SAVE  FI 
FETCH  C 
DEFINE 

ADO  OTH 
FNO  OF 
ITERATE 
GO  ON  I 
SET  UK 
QUIT  IF 
ZERO  CO 
START  w 

7600000 
RESULTS 
ROUTING 
SAVE  UT 
SAVE  FI 
FETCH  C 
DEFINE 

ADD  OTH 
END  OF 
ITERATE 
GO  ON  I 
XSET  UP 
QUIT  IF 
ZERO  CO 
START  « 

0000000 
RESULTS 
ROUTE  E 
SAVE  IN 
ENABLE 
INTERCH 
ANO  REG 


RESUL 
ITERAT 

ONCE 
F  NOT 

EXPQNE 
SUE* 

UNT. 

ITH  FI 

181 
FROM 

OEFIN 
HER  OP 
RSTOAP 
ONSTAN 
(RGS*S 

ER  OPE 
ITERAT 

ONCE 
F  NOT 
EXPQNE 

SUE- 
UNT. 
ITH  Fl 

i  a  i 

FROM 
OEFIN 
HER  OP 
RSTOAP 
ONSTAN 
(RGS*S 
ER  OPE 
ITERAT 

ONCE 
F  NOT 
EXPON 
SUE* 
UNT. 
ITH  FI 
181 

FROM 
NO  ARO 

REGIS 
PES  0 

ANGE  R 
ISTER 


T, 

ION. 

FOR  EACH  RON. 

DONE. 

NT  FOR  INVERSE. 

a. 

RST  ROw  OF  DATA. 

PREVIOUS  STEP. 

E. 

ERANO, 

ERANO. 

TS  (W  S). 

IN). 

RANO* 

ION, 

PER  ROW, 

DONE. 

NT  FOR  INVERSE. 

16, 

RST  R0W  OF  DATA, 

PREVIOUS  STEP. 

E 

ERAND. 

ERANO. 

TS  (*  S). 

IN), 

RANO. 

ION. 

PER  ROW. 

DONE. 

ENT  FOR  INVERSE. 

32. 

RST  ROW  OF  DATA, 

PREVIOUS  STEP. 

UNO  32, 

TER  S, 

THRU  3l, 

EGISTER  % 

A. 


-1+0- 


NE64I 
F128| 


NE128l 

P256I 

F256j 


SKIP1  I 

NE256» 

P512I 
F512, 


SKIP2J 


NE512I 

P1024J 
F1024I 


SETE 

SETE1 

STS 

LDS 

LDA 

COMPL 

AON 

STA 

TXLTM 

CTS8F 
LIT(2 

Jump 

sutc 

caoo( 

LDS 

MULTA 

AUTC 

TXLTM 

CTSBF 

LIT(2 

JUMP 

SUTC 

CADO( 

RESET 

LOS 

MULTA 

ALITf 

SKIPT 

CHANG 

SETFF 

SKIP 

ALIT( 

TXLTM 

CTSBF 

LIT(2 

JUMP 

SLIT( 

CADO( 

LIT(L 

RESET 

LOS 

MULTA 

ALIT( 

SKIPT 

CHANG 

SETFF 

SKIP 

ALIT( 

TXLTM 

ALIT( 

ALIT( 

TXLTM 

CTSBF 

LIT(2 

JUMP 

SLIT( 

CADDC 

LIT(L 

KESET 

LOS 

MULTA 


EXMULT 


(COUNT) 
(SIZE) 

) 

COUNT) 
COUNT) 

QDSUR 
COUNT) 

(COUNT) 

(SIZE) 

) 

COUNT) 
COUNT) 
FF 

ODSUB 
COUNT) 

EQUAO 


COUNT) 
(COUNT) 
(SIZE) 
) 

COUNT) 
COUNT) 

OOP) 
FF 

ODSUB 
COUNT) 

EQUAD 


COUNT) 

(LOOP) 

LOOP) 

COUNT) 

(COUNT  ) 

(SIZE) 

) 

CUUNT) 
COUNT) 
OOP) 
FF 

ODSUR 


E.OR. 

E'OR* 

0<COU 

SAJ 

SIN*3 

SIN+3 

0(COU 
0(COU 
.F64» 
33.ME 
PAIR( 

nuT> 

■  01 
SSIZE 

SlN  +  4 

1  » 
.11 
•  Fl?R 
32.  NE 

PAIR( 
OUT  l 
.0; 
SSTZE 

J 
SIN  +  5 

2; 

■  11 

.SKIP 
SIN  +  5 

I 
»F256 

■  II 

.P256 
31»nE 
PAIR( 
UUTl 

■  o; 

SSIZE 
1*1.0 

; 

6  +  SIN 
4| 
-2) 
.SKIP 
6  +  SIN 

1 
.F51? 
a-3J 

»PS12 
-2; 

*5l 

»PSl2 
30.  Nt 

PAIr?( 
OUT; 

■  01 
SSIZE 
1.3,0 

i 
8  ♦  S  I  N 
HI 


-El 
"El 
NT)I 

I 
I 

NT)! 
NT)I 

64, 
32.U»3?.0 


I 

uai 

64,0*64.0 


U 

I 

I 
I 

2b6; 

12«.u»i2* 


i 

(LOOP)  l 


2» 

(LOOP)  I 

I 
I 


I 

512: 

2bft. J,?56 


; 
I 

(LUUP)I 


IENA 
X 

»SAV 
*SAV 
XFET 
IDEF 
XADD 
*END 
«ITE 
XfiO 
)IXSE 
XQui 
XZER 
XSTA 

XFET 
XOEF 
*GU 
«ITE 
XQU 
)  JXSE 
XQUI 
XZER 
XSTA 
XoEF 
XFET 

xnEr 

*INC 
»CHE 
*OEF 
XDEF 
IINN 
XINC 
XlTE 
XGO 
.0)1* 
*QUI 
*ZER 
*STA 

*l_00 
XoEF 
XFET 

xoEf 

*INC 
*CHE 
*DEF 

xoEf 

*INN 

*INC 

<G0 

XRES 

XSET 

*r,o 
xr,u 

.0)11 

xqui 

XZER 
XSTA 
XLUO 
XDEF 
XFET 
XDEF 


BLE  ALL  PES. 


E  UTHER 
E  FIRST 
CH  CONS 
INE  (RG 

OTHER 
OF  ITE 
RATE  ON 

ON  IF  n 
T  UP  EX 
T  IF  SI 
0  COUNT 
RT  WITH 

CH  CONS 
INE  (D 
BY  PAIR 
RATE  HN 
ON  IF  N 
T  UP  EX 
T  IF  SI 

n  count 

RT  WITH 
INE  (TF 
CH  CONS 
INE  (oA 

REMENT 
CK  INNE 
iNt-MOD 
INE  (TF 
ER  LOOP 
RFMENT 
RATE  ON 
ON  IF  N 
SFT  UP 
T  IF  SI 
0  COUNT 
RT  WITH 

P  COUNT 
INE  (TF 
CH  CONS 
INE  (DA 

RFMENT 
CK  INNE 
InE-MOO 
INE  (TF 
ER  LOUP 
REMENT 
BACK  FQ 
ET  LOOP 

COUNT 
RACK  FO 
ON  IF  N 
SET  UP 

T  IF  SI 

0  COUNT 
RT  WITH 
P  COUNT 
iNt  (TF 
CH  CONS 
INE  (OA 


OPERAND. 
UAPERAND. 
TANTS  (W  S). 
S#SIN), 

OPERAND* 

RATION. 

CE  PER  ROW, 

OT  DONE. 

PONENT  FOR  INVERSE. 

*E-64. 

FIRST  RUW  OF  DATA. 
TANTS  (W  S). 
ATA  +  1  )*SlN  +ANO-  DATA 
S  OF  ROWS, 

CE  PEp  PAIR  Of  ROWS* 
OT  DONE. 

PONENT  FOP  INVERSE. 
ZE-12B. 
. 

FIRST  RUw  OF  DATA. 
FF«0)  LOOP  CONTROL. 
TANTS  (*  S). 
TA*2)*slN  «.ANn-  DATA, 

COUNT. 

R  LOOP. 

IFY  CONSTANTS. 

FF«1)  LEAVE  NEXT  TIME 

. 

COUNT. 

CE  FOR  EACH  4  ROWS. 

OT  DONE. 

EXPONENT  FOR  INVERSE. 

ZE-256. 

VlRST    ROw    OF    DATA. 

ER. 

FF«0)  LOOP  CONTROL. 

TANTS  (W  S). 

TA*4)*SIN  +AND-  DATA. 

COUNT. 

R  LOOP, 

IFY  CONSTANTS. 

FF-1)  LEAVE  NEXT  TIME 

. 

COUNT. 

R  NEXT  PAIR  OF  ROWS. 

COUNTER. 
FOR  NEXT  OATA  SET. 
R  NEXT  DATA  SET. 
OT  DONE. 

EXPONENT  FOR  INVERSE, 
ZE-512. 
. 

FIRST  ROW  OF  DATA. 
ER. 

FF«0)  LOOP  CONTROL. 
TANTS  (W  s). 
TA+8)*SIN  +ANn-  DATA. 


-1+1- 


S*IP3« 


NE1024» 

P204BI 
F2048i 


SKIP4I 


NE2048' 
Pa096i 


SKIP5I 


OuTl  I 


A|_IT(C 

SKIPT 

CHANGE 

setff 
skip 

ALIT<C 
T*LTM( 
*LIT(L 
ALIT(C 

TX|_TM< 

ctsbfc 

LIT(2) 

JUMP 

SLITCC 

CADOCC 

LIT(LU 

HF.SETF 

LOS 

MULTAU 

ALITCC 

SKIPT 

CHANGE 

SETFF 

SKIP 

ALITCC 

TXLTM< 

ALIT(L 

ALITCC 

TXLTMC 

CTSBFC 
LIT(2) 

JUMP 
SLITCC 
CADOCC 
LIT(LO 

RESETF 

LOS 

MULTAD 

AlIT(3 

SkIPT 

CHANGE 

SETFF 

SKIP 

ALITCC 

TXLTM< 

ALIT(L 

ALITCC 

TXLTMC 

LIT(2) 

CTSBFC 


OUNT) 

Quad 


HUNT  ) 

LOOP) 
OOp) 

OUNT) 

COUNT) 

SIZE) 


OUNT) 
OUNT) 
OP) 
F 

OSUR 
OUNT) 

QUAD 


OUNT) 
LOOP) 
OOP) 

OUNT) 

COUNT) 

SIZE) 


OUNT) 
OUNT) 
OP) 
F 

dsur 
) 

QUAO 


OUNT) 

LOOP) 

OOP) 

OUNT) 

COUNT) 

SIZE) 


•  41 
.SKIP 
8*SIN 

I 

•  F102 
—  7» 

,Pl02 
.-«! 

-in 

#P1()2 
29»nE 
pAI^( 

OUTi 
■  01 
SSIZE 

1  p  7  ,  0 


3» 

( LOUP) I 

4) 
41 


INCREMENT  COUNT. 
*CHECK  INNER  LOOP, 

IdEFINE-NODIFV  CONSTANTS, 

tntriwi  (Trrr«l)  liavc  NtxT  time 


1024» 
512.0.5 


41 
NCLOOP) 

i8| 

81 


.8» 

:2U4rt» 

1U24.Q. 


12+SINCL00P)! 
161 
ubl 
.SKIP 
12*SI 

I 
•  F204 

■  -IS) 
»P204 

■  •8| 

■  23| 
.P204 
28.NE 
PAI^( 
(JUT, 
•01 

SSI  ZE  i 
1*15.0} 

. 
20*SINClOUP)» 
32l 

■  16} 

♦SKIP5J 
204-SINCLnUP) 

I 
»F4o96i 

■  -31  J 
.P4096J 
«-l6l 

■  47, 
.P4096J 
PAIR(2U48.0. 
15.RETURN1 


1INNER  L00P» 

•INCREMENT  COUNT, 

XGQ  RACK  FOR  NEXT  PAjR  OF  HOWS. 

fPtSET  LOO?  cOuNTEH. 

XSET  COUNT  TOR  NEXT  DAT*  SET. 
*G0  RACK  FOR  NEXT  OATA  SET. 
«G0  ON  IF  NOT  OONE. 

!2.0)hset  up  Exponent  for  inverse, 

IQUIT  IF  DONE  CSIZt»l024), 

XZERO  COUNT, 

XSTART  WITH  FIRST  ROW  0^  DATA. 

XLOOP  COUNTER, 

DINNER  LOOP  CONTROL,  CTFFF"0>, 

XFETCH  CONSTANTS  (N  $>. 

XDEFINE  (DATA*16)*SIN  *AND-  DAT* 

^INCREMENT  COUNT, 

ICHECK  INNER  LOOP, 
»   SQEFINE-MODIFY  CONSTANTS, 

XnEFlNt  (TFFF.l)  LEAvE  NEXT  TIME 

tlNNER  LOOP. 

INCREMENT  COUNT, 

XGQ  RACK  FOR  NEXT  PAIR  OF  ROWS. 

XRESET  LOOP  COUNTER, 

XSET  COUNT  FOR  NEXT  DATA  SET, 

*G0  RACK  FOR  NEXT  DATA  SET. 

*GU  ON  IF  NOT  DONE. 
1024,0)1      XSET  UP  EXPONENT. 

XQUIT  IF  DONE  <SIZE»?0a8). 

XZERO  COUNT, 

*START  WITH  FIRST  HOW  nF  DATA, 

XLOOP  COUNTER, 

XDEFINE  (TFFF-O)  LOOP  CONTROL. 

XFETCH  CONSTANTS  CW  S), 

XDEFINE  (DATA+32)#SIn  *AND-  DATA 

XINCRFMENT  COUNT, 

ICHECK  INNER  LOOP, 
I   *DEFlNE-MnOIFY  CONSTANTS. 

XnEFlNE  (TFFF«i)  LEAVE  NEXT  TIME 

IINNER  LOOP. 

*INCRFMENT  COUNT, 

XGO  RACK  FOR  NEXT  PAIR  OF  ROWS. 

XRESET  LOOP  COUNTER. 

XSET  COUNT  FOR  NEXT  DATA  SET. 

XgO  RACK  FOR  NEXT  DATA  SET.. 
2048,0)1      XSET  UP  EXPONENT. 

XI.EAVE  IF  NOT  INVERSE  pFT, 


THE  FOLLOWING  COUE  IS  EXECUTED  IF  AN  INVERSE  TRANSFORM  IS 
REQUESTED.   THE  ONLY  OPERATION  REQUIRED  IS  THAT  EACH  ELEMENT  OF 
THF  SOLUTION  M 1 1 s T  BE  UlVIUEn  BY  N(SIZE).  THE  NUMBER  OF  DATA  POINTS 
IN  TME  FUNCTION,   SINCE  N  Is  A  ROWER  OF  TWO  THIS  CAN  BE  DONE  BY 

SUBTRACTING  LOG(RASE  2)^i    FHQH  EACH  FXPONFNT.   THE  PROPER  EXPONENTS 
AHr  CREATED  AT  A  1  I  HE  WHEN  NCSIZE)  IS  KNOW  EXACTLY  AND  THE*  ARE 
PUT  IN  ACAR2CLHHP)  JUST  HEFORt  EXECUTION  IS  TRANSFEREO  TO  OUT. 


SLITCCOUNT)     aQJ 
CADDCCOUNT)     SIZE» 


XZERO  COUNT. 
XSTART  WITH  FIRST 


ROW  OF  DATA. 


-1+2- 


INVERSE 

LOA 

O(COUNT) J 

SBEX 

SC2I 

TXLTM(COUNT) 

.INVERSE! 

RETURN! 

LDL(3) 

*02| 

LOL(O) 

*D3l 

LOKl) 

$04> 

EXCHL(3) 

SlCRj 

ENO 

FFT1. 

XFETCH  DATA. 

XDIVIDE  BY  SIZE  (NO  uNnERFLOW). 

XDO  FOR  ALL  ROWS  OF  DATA. 

*FETCH  RETURN  AqdRESs. 

IRESTORE  ACARO. 

^RESTORE  ACAR1, 

XRETURN  TO  CALLING  PROGRAM. 


43- 


812) 


S  pESt-T  LIST! 
BEr,lN 

FILL       l28t 
SIwUnTRYHiDATA  ( PA  IRC  1.0* 1, 0 )• PA IR(, 707 10678 12..707 1067 
PAIR(  1.0* -1.0). PA  I R<-. 70710678 12*. 707 1067812)* 
PAIR(-1.0»-1.0).PAIR( -.707 10678 12. -.707 1067812), 
PA  I R( -1.0»l,n), P A  IRC. 7071067812.-. 707 10 A78t 2 ))8. 
(  ^A I R( 1, 0» 0.0), PA IR<« 9238795 325 », 38268  3 A  324). 

PMR(  0.707106T8i?»  0.  7071  0678  1  2) .  PA  IRC  0.  3826834324,  0,9238795325)* 
PMR(0»0»  1. 0) » P At R<"« 38268 34 324*.92 38795325)* 
PAIR( -0,70  7106781 2*  0. 707106781 2), PA  I R( -0,9238795325, 
PAlR(-l,0»0.0),PAlR(".9238795  3?5»-*38268  34  324). 
PA  I  R^ -.  7  o7  t06  78l  2, -.7  07  10678  12).  PA  IRC  3826834  324* -.9238 

PAIR<0.0»-1«0).P A TH(» 3826834 324»-,923875 325). 
PaIPI  0,707l0678l?»-0.707l0678l2)»PAIR(  0 . 923879532s, -0 , 
( KA I R( 1. 0» 0«0). pa  I R(» 9807852804* .19509  0322)* 


0.382683432M* 

795325). 

3826834124))4. 


PAIR(  0,923879s??5»  0.38?68343?4 
PaIr(  0.70710678!?.  0.7o7lo678i2 
PAIR(  0,38268343?4»  0. 92387953?5 

KAIhfO.O*  1.0),  PAIRC  195090322* 
PAlR(-0.38?683/l3?a»  0.92387953?5 
PAlR(-0.7o7lo6T8i?»  0.7o710678i2 
PAlR(-0,923879s3?5»  0,38268343?4 

KAlRf-l.O,0.0).PAlR(-. 980785 280 
PAlR(-0, 923879s? ?5»"0. 38 268343?4 
PaIr(-0.7o710678i2'-0.7o710678i2 
PAlR(-0,3a2  683«3?A»-0,92  387953?5 

HAIR(0,0»-1.0)»PAlRC 195  09  0322* 

PAIr(  0.38 2683a 3 ?4#-0. 92 38 795325 

PaIrC  0,7o7lo^r8i2*-o.7o7lo678i2 

PaIR^  0,<?23  87953?s»-n»382  6834  324),p 

^AIR(1,0»0.0),PAIR(*9951847267» 
PAIr(  0,980^85?804»  0.1950903220 
PAIe,(  0.92387953?5»  0.38?68343?4 
PAIR(  0.83l469^1?3,  0,5555702330 
PAIR(  0,707l067fli?»  0.7071067812 
PAIr(  0.55557o?3io»  0. 8314696123 
PaIr(  0,3fl2683/j3?a»  0.92387953?5 
PAIR(  0, l95090^2?0»  0,9807852804 

>*AI«f0,0»1.0).PATR(".  0980171403 
PAIrC-O. 1950903220*  0.98o78528o4 
PaIr(-0.38  26B3«3?4»  0.9 2 38 795 3?5 
PAIR(«0,555570?3?0,  O.B^H^l?3 
PAIR(-0.707106781?.  0,70710678i2 
pAI»(-0.83l46961?3»  u. 5555702330 
PairC-0.92  38  79s3?s»  0.38 268 34 3 ?4 
PMR( -0,980  785?  804.  0,1 9  50  903220 

KAIR("1.0*0.0)*PAIR(-*9951847  26 

PAIR(-O.9  8o78  5p8o4»"0. 195  09  03220 
pAlR(- 0.9 2 3»79S3?5 '"0.38  268 34 324 

PA  I  R( -0,831  a  69A|i  ?3.-0,  5555  70233  0 
PAIRC- 0,7071 0^78 12* -0.707 10678 12 
PA  I r( -  0. 55557 o?330»" 0.83 14696123 


)»PAIR<  0.8314696123 
).PAIR(  0.5555702330 
).PAIR{  0,1950903220 
.9807852804)* 

)»PAlR(-0. 5555702330 
).PAlR(-0. 8314696123 
),PAlR(-0, 9807852804 
4.-. 195090322). 
).PAIR("0. 8314696123 
)*PAlR(-0. 5555702330 
).PAlR(-0. 1950903220 
-.9807852804). 
)»PAIR<  0.5555702330 
)*PAIr(  0.8314696123 
AIR(  0.980785?804,-0 
.0980171403)» 
).PAIR(  0.9569403357 
).PAIR(  0.8819212644 
).PAIR(  0,7730l04534 
).PATR(  0.6343932842 
)iPAlR(  0.4713967368 
)»PAIR<  0.29028467^3 
).PAIR(  0,0980171403 
..9951847267), 
)*  PAlR<-0, 2902846773 
)»PaIr<-0. 4713967368 

).PAlR(-0, 63439328*2 
).PAIR("0. 7730104534 
).pAlR(-0. 8819212644 
).PAIr(-0. 9569403357 

)»PAlR(-0, 995184 7267 
7.-. 0980171404  ). 
)»PAIR("0. 9569403357 
).PAIR("0. 8819212644 
),PAlR(-0, 7730104^34 
).  PAIRCO.  6343932842 
)»PAlP(-0. 4713967368 
)»PAlR<-0. 2902846773 


P  AIr(-0.  38  ?68  343->4»- 0.9  238795  325 
PA IR^-.l 950903220,-, 080^852804 ),PAIR( -.09801 7l*03.-.995l 
KAlR(0.0*-1.0).PAlRC.0980l7l403»-# 9951847267). 
pair(  0, I950903??0»"0.98o78528o4), PAIR(  0,290284677  3, 


0.5555702330) 
0.8314696123) 
0.9807R52804) 

0.8314696123) 
0.55557n2330) 
0,l9S0903220) 

•0.5555702330) 
•0.8314696123) 
-0,9807A52804) 

-0.8314696123) 
-0.5555702330) 
1V50903220))2, 

0.29o2846773) 
0.4713967368) 
0,6343932842) 
0.7730104534) 

0.88l9?i2644) 
0.9569403357) 

0.995M47267) 

0.9569403357) 
0.88l92l2644) 
0,7730104534) 
0.6343932842) 
0.4713967368) 
0.2902846773) 
0,0980171403) 

-0.2902846773) 
-0. 4713967368) 
-0,6343932842) 
-0.7730104534) 
-0.8819212644) 
-0.9569403357) 

847267)» 


•0.9569403357). 


00000090 
00000100 
00000110 
00000200 
00000300 
00000400 
00000500 
00000600 
00000700 
00000800 
00000900 
00001000 
0000110O 
00001200 
00001300 
00001400 
•00001500 

•  00001600 

.00001700 

OOOOlROO 
•00001900 
•00002000 
.00002100 

00002200 
•00002300 
»  00002400 
•00002^00 

00002600 

•  00002700 

•  00002*00 
00002900 
00003000 

•  00003100 

•  00003200 
.00003300 

•  00003400 
• 00003SOO 

•  00003600 
»00003r00 

00003800 
»  00003900 

•  00004000 
#00004100 

.00004200 
»  00004300 

•  00004400 
.00004^00 

00004*00 
• 0000470O 
• 00004A0O 
.00004Q00 

•  00005000 

•  00005100 

•  00005200 
0000s300 
00005400 
00005500 
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